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INTRODUCTION 


Since  the  July,  1980  conception  of  the  Kansas  Wind  Resource 
Assessment  Project  at  Kansas  State  University,  wind  data  at  eight 
Kansas  locations  has  been  collected  using  stand-alone  wind  data 
acquisition  systems.  After  data  was  collected  and  stored  on  tape 
in  the  field,  it  was  transferred  to  KSU  by  mail,  processed,  and 
stored  on  the  Wind  Project  Data  General  Nova  computer.  Now  that 
a  substantial  data  base  has  been  established,  studies  of  wind 
energy  production  can  now  be  made.  This  thesis  incorporates 
three  years  of  wind  and  utility  data  from  1980  to  1982,  and 
analyzes  the  effects  that  wind  power  generated  during  this  time 
could   have   had   on   conventional   utility  operation. 

To  analyze  the  effects  of  injected  wind  power,  five  major 
calculations  are  used.  In  each  calculation,  ramping  is  a  term 
applied  to  the  rate  of  change  of  power  output,  and  is  the  princi- 
pal figure  of  merit  used  throughout  this  thesis.  The  major 
calculations  consist  of  a  24  hour  average  and  worst  case  ramping 
day,  a  12  month  average  and  worst  case  ramping  year,  an  average 
and  worst  case  ramping  weekend  and  weekday,  a  distribution  of 
ramping  rates,  and  a  percent  energy  loss  curve  due  to  upwards 
ramping  control.  Each  of  these  calculations  is  discussed  in 
detail  below  within  Chapter  III,  and  aid  in  detecting  the  effects 
that    injected  wind  power  has   on   the   utility. 

In  addition  to  the  induced  ramping  effects  on  utility  opera- 
tion,  a  correlation  study  is  also  given.     This  brief   study  uti- 
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lizes  autocorrelation  and  crosscorrelation  functions,  and  allows 
detection  of  wind  power  and  Kansas  load  demand  patterning. 

Next,  to  establish  a  base  from  which  to  work,  definitions 
for  analysis  are  given  to  acquaint  the  reader  with  three  terms 
used  throughout  this  thesis. 

Percent  Penetration  -  the  fraction  of  rated  wind  power  out- 
put to  rated  utility  power  output  (assumed  to  be  3000  MW). 

Farm  -  a  single  wind  site,  possibly  consisting  of  numerous 
wind  turbines,  all  under  identical  wind  conditions.  Since  all 
turbines  at  a  site  are  assumed  under  identical  wind  conditions, 
the  terms  wind  farm  and  wind  turbine  will  be  interchangeable. 

Wind  Array  -  two  or  more  farms  covering  a  state  wide  area. 

As  noted  above,  a  3000  MW  penetration  base  is  used.  This  is 
done  since  the  plot  of  utility  power  production  appears  to  peak 
around  3000  MW  (a  good  round  figure)  for  Kansas  Power  and  Light 
and  Kansas  Gas  and  Electric  combined  power  production.  Also, 
these  utilities  were  used  since  they  represent  approximately  75% 
of  Kansas  power  production,  and  had  demand  data  that  was  conven- 
iently available. 

Next,  to  establish  a  basic  knowledge  of  the  four  sites  used 
in  this  study,  a  brief  description  for  each  site  is  given,  while 
approximate  geographical  locations  of  each  site  are  shown  in 
Figure  1.1. 

Tuttle  Creek,  the  first  site  established  under  the  Wind 
Resource  Assessment  Project  at  KSU,  is  located  on  the  Corps  of 


Engineers  tower  at  Tuttle  Creek  reservoir  near  Manhattan,  Kansas. 
It  was  established  July  3,  1980  and  has  a  base  elevation  of  1250 
feet  above  sea  level.  The  tower  lies  to  the  west  end  of  Tuttle 
Creek  dam  on  a  well  exposed  hill,  while  normal  water  level  of  the 
reservoir    is    1075    feet    above    sea    level. 

Wright,  which  only  has  a  maximum  30  m  tower  height,  is 
located  on  the  KTVC-TV  microwave  relay  tower  at  Wright,  Kansas. 
The  site  was  installed  August  15,  1980,  and  has  a  base  elevation 
of  2530  feet  above  sea  level.  This  tower  is  located  at  the 
south-east  corner  of  Wright  near  a  school,  is  open  to  the  south 
and  east,    and  has   a   grain  elevator   1/3   mile    to   the  northwest. 

Plainville.  located  on  the  Dale  Roll  farm  2  miles  west  of 
Plainville,  Kansas,  is  positioned  on  an  Oil  Field  Communication 
tower,  and  was  installed  October  11,  1980.  It  has  a  base  eleva- 
tion of  2160  feet  above  sea  level,  and  is  well  exposed  from  all 
directions  except  for  a  house  and  trees  100  yards  east  of  the 
tower. 

Finally,  Atlanta,  established  September  25,  1981,  is  located 
1  and  3/4  miles  north  of  Atlanta.  The  tower  has  a  base  elevation 
of  1430  feet  above  sea  level  and  lies  on  a  hill  which  is  well 
exposed   from    all    directions. 

To  conclude  the  introduction,  special  note  is  made  regarding 
the  plots  used  in  subsequent  analyses.  Since  over  120  plots  were 
generated  for  this  study,  only  samples  from  the  original  number 
are  used  as  examples  within  this  thesis.  If  additional  plots  are 
needed  to  verify  any  results  below,  refer  to  the  Kansas  Wind 
Resource  Assessment  Project   at   Kansas    State   University. 


CHAPTER  II 
LITERATURE  REVIEW 

Much  research  has  already  been  done  in  the  area  of  wind 
generation  effects  on  utility  systems.  In  fact,  most  studies  to 
this  date  analyze  the  effects  that  thunderstorm  induced  wind 
generation  have  on  utility  operation.  The  major  analyzing  method 
used  in  these  studies  involves  modeling  wind  behavior  and  utility 
operation,  with  some  reports  using  actual  measured  wind  speeds  to 
simulate  wind  turbine  power  output.  Finally,  the  reports  then 
study  the  effect  that  injected  wind  generation  has  on  area  con- 
trol   error    or   ACE. 

In  one  study  by  Schlueter,  Park,  Modir,  Dorsey,  and 
Lotfalian  [8],  worst  case  effects  of  wind  farm  power  production 
due  to  thunderstorm  fronts  were  analyzed.  The  paper  utilized  a 
modeled  wind  speed  profile,  the  Mitchell  model  [10],  to  analyze 
how  this  dramatic  change  in  power  affected  the  overall  area 
control  error.  One  of  the  major  points  of  the  paper  stated  that 
given  a  high  enough  wind  power  penetration,  the  rapid  increase 
and  decrease  in  wind  power  production  due  to  thunderstorms  could 
cause  a  violation  of  utility  performance  requirements,  an  example 
being  ACE  equaling  zero  at  least  one  time  in  all  ten  minute 
periods.  However,  the  paper  states  that  this  problem  can  be 
controlled  by  limiting  the  total  wind  generation  capacity  to 
less  than  the  automatic  load  following  capability  of  automatic 
generation  control,  limiting  echelon  [8]  penetration  between  2 
and   3%,    and  by  shutting   down  portions   of   a  farm  before  a   thunder- 
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storm    pass. 

In  a  second  paper  [9],  Schlueter,  Park,  Lotfalian, 
Shayanfar,  and  Dorsey  expanded  on  the  paper  given  above.  Since 
a  definite  problem  in  load  following  capability  could  exist  under 
adverse  conditions,  the  authors  investigated  modifications  to 
power  systems  operation  to  compensate  for  significant  wind  power 
penetration.  These  modifications  consisted  of  setting  higher 
requirements  for  spinning  reserve,  load  following  capability,  and 
unloadable  generation.  Thus,  the  problem  of  increased  ACE  could 
be  reduced  by  modifying  parameters  on  existing  generation  meth- 
ods, and  would  be  desirable  since  implementation  of  a  complex 
control  system  for  the  utility  control  of  dispersed  wind  genera- 
tion would  be  costly.  However,  when  such  standards  are  set 
higher,  production,  operation,  and  maintenance  costs  would  like- 
wise increase  since  units  not  otherwise  required  would  be  run- 
ning. 

In  their  most  recent  paper  [10],  Schlueter,  Park,  Lotfalian, 
Shayanfar,  and  Dorsey  introduce  three  more  methods  for  reducing 
thunderstorm  induced  power  changes  from  large  wind  tubine  arrays. 
These  methods  include  selection  of  a  wind  turbine  model  for 
particular  sites,  selection  of  an  appropriate  siting  configura- 
tion, and  wind  array  controls.  If  these  methods  are  used,  then 
previous  wind  farm  requirements  set  in  [8]  would  not  depend  on 
the  capacity  of  the  farm,  but  on  the  wind  turbine  selection  and 
siting  configuration,  while  coordinated  blade  pitch  controls 
would  be   used   to  reduce    the    effective    farm   penetration   level. 


Also  of  interest  is  the  effect  that  small  wind  turbines  have 
on  utility  generation  performance.  In  a  paper  by  Curtice  and 
Eeddoch  [11],  a  simulation  model  of  a  utility's  automatic  genera- 
tion control  was  modified  by  including  synthesized  data  repre- 
senting an  aggregate  output  of  small  turbines.  Since  small  wind 
turbines  have  less  inertia  than  large  turbines,  their  output 
power  variations  would  tend  to  be  more  frequent.  In  essence,  for 
turbines  with  variations  greater  than  0.01  cycles  per  second,  an 
added  noise  component  of  relatively  higher  frequency  would  be 
injected  into   the  power  system.  This   component    is   likely  to  be 

uncontrollable,  and  if  large  enough,  could  cause  the  quality  of 
the  utility's  system  performance  to  decrease.  That  is,  higher 
ACE  values,  longer  times  between  zero  crossings,  etc.  would 
result. 

In  summary,  the  papers  discussed  above  show  how,  given 
significant  penetration,  wind  generation  can  produce  undesirable 
effects  on  a  utility's  load  following  capability.  Also,  several 
methods  are  mentioned  by  which  rapid  changes  in  power  output 
from  large  wind  turbine  arrays  could  be  reduced,  and  how  modifi- 
cation of  power  system  operation  could  be  used  when  there  is 
significant  wind  power  penetration.  In  this  thesis,  the  two 
methods  used  to  aid  in  controlling  negative  effects  due  to  injec- 
ted wind  generation  will  be  upwards  ramping  and  wind  array 
controls. 


CHAPTER   III 

METHODS,    PROCEDURES,    AND  DATA 

Within  this  chapter,  defining  relations  and  data  validity 
are  discussed  to  give  insight  into  the  analysis  of  wind  power 
effects  on  conventional  generation  for  this  thesis.  The  main 
areas  of  discussion  involve  calculations  for  wind  power,  ramping 
rates,  ramping  distributions,  ramping  control,  and  correlation, 
with  an  additional  section  on  the  validity  of  the  wind  and  utili- 
ty  sample    data    sets. 

Wind  Power  Production 

In  calculating  wind  power  from  mean  hourly  speeds,    the  model 
used    will    be    that    in   Figure    3.1.      This    model    incorporates    a 
linear  relation  between  wind  speed  and  power  output   from    the   cut- 
in   wind    speed    (u    )    to    the    rated   wind    speed    (u    ).   The   relation 
c  r 

between  power  and  wind  speed  is  actually  a  cubic  relationship, 
but  as  stated  in  [5],  low  efficiencies  at  low  wind  speeds  help  to 
linearize  the  power  output  curve.  This  is  also  verified  by  the 
actual  power  output  operation  of  a  MOD-2  wind  turbine  shown  in 
Figure  3.2.  This  thesis  will  base  power  production  calculations 
to   simulate    turbines    similar    to  the   MOD-2. 

Also,     considering    the    furling   wind    speed    (u_),    no   action    is 
taken    to    shut    down    power    production   during    calculations    for    two 
reasons.      First,    the    turbine  model   used   in   this    thesis    incor- 
porates u     values  of  5,  6,  and  7  mps  and  u     values  of  10,  12,  and 
c  r 

14   mps    respectively.      Since    these   values    closely   follow    those    of 
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Figure  3.1.  Wind  turbine  power  output  vs.  windspeed. 
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Figure  3.2.  Performance  curve  for  MOD-2  wind  turbine 
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the  MOD-2,  which  can  have  a  60  mph  furling  wind  speed  [4],  and  no 
wind  data  used  in  this  report  exceeds  50  mph,  no  furling  condi- 
tion is  ever  reached.  Second,  the  wind  speeds  used  in  calcula- 
ting modeled  wind  power  production  are  mean  hourly  speeds.  The 
collection  time  period  for  this  data  is  not  short  enough  to 
detect  gusting  wind  speeds  in  excess  of  60  mph,  and  thus  need  for 
such  detection  in  power  calculations  would  be  unnecessary.  The 
above  two  arguments  do  not  neglect  the  fact  that  a  wind  speed 
above  60  mph  could  have  occurred  during  valid  data  collection, 
but  suggest  that  detection  of  such  conditions  could  only  occur  if 
such  wind  speeds  persisted  over  an  hourly  period  of  data  collec- 
tion. 

For  the  model  shown  in  Figure  3.1,  a  more  precise  expression 
for  the  average  power  output  of  the  turbine  is 

u  <  u 
c 

uc  i  u  i.  ur        (3.1) 

u  >  u 

r 

where   P      is    the   electrical   power   out,    P„   is    the   rated   power,     and 
e  K 

u  is  the  mean  hourly  wind  speed.  A  sample  plot  of  power  output 
using  this  equation  can  be  seen  in  Figure  3.3.  This  output  power 
landscape  is  calculated  from  wind  speeds  measured  at  the  Tuttle 
Creek  50  m  tower.  Two  successive  blow-ups  of  this  landscape  show 
the  hourly  power  output  over  a  120  hour  period  which  corresponds 
to  April  6,    1982  at  8  am   to  April  11,   1982  at  12  pm.      Note    that 
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the   scale  on  the  horizontal   axis  of  Figure  3.3    is   in  hours,   and 
corresponds    to   the   hour    in   year    1982. 

Another  important  note  to  make  concerns  the  injection  of 
wind  power  into  the  utility  grid.  This  will  be  done  by  first 
assuming  that  the  utility  power  data  matches  Kansas  load  demand. 
In  using  this  assumption,  since  coincident  wind  and  utility  power 
production  will  be  adding  together  to  meet  Kansas  load  demand,  to 
find  the  resultant  utility  power  output,  wind  power  will  be 
subtracted  from  the  Kansas  load  demand.  An  example  plot  of  this 
operation  can  be  seen  in  Figure  3.4.  Note  that  this  figure  has 
the    coincident    wind    power    production    in    the    lower    plot. 

Ramping  Rates 

Calculating  the  ramping  rates  from  the  wind  and  utility 
demand  data  will  be  done  by  utilizing  their  power  output  as  the 
average  power  over  an  entire  hour.  Thus,  if  KPL  and  KGE  satis- 
fied a  Kansas  load  demand  of  1000  MWhrs  for  a  particular  hour  of 
the  day,  the  average  power  output  would  be  1000  MW  over  the 
entire  hour.  So,  on  an  hour  by  hour  basis  using  average  hourly 
power   output,    ramping  rates  are  calculated  using 


Pave   -  Pave 

Ramping  Rate   =  l     MW/min  (3.2) 

60 


where  Pave^^  and  Pave2  are  consecutive  average  hourly  powers,  and 
60  has  the  dimension  of  minutes.  For  increasingly  smaller  time 
intervals,  this  expression  would  approach  the  derivative  or  in- 
stantaneous rate  of  change  of  the  electrical  power  output.   In  a 
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more  general  form. 


P.  ,  -  P. 

R.  =  -* L     MW/min  (3.3) 

J      60 


where  R  is   the  average  hourly  ramping  rate,   and  P  is  the  hourly 
average    power   production. 
24  Hour  Ramping   Calculations 

24  hour  ramping  rates  are  best  described  as  24  hours  in 
which  three  years  of  data  have  been  condensed.  That  is,  each 
hour  for  the  entire  period  is  summed  or  compared  with  its  fellow 
hour  of  another  day.  The  final  results  from  the  24  hour  calcu- 
lations give  mean  ramping  rates  for  a  24  hour  day,  standard 
deviations  of  the  periods  ramps  with  their  corresponding  mean 
ramp,  and  maximum  and  minimum  ramps  for  a  worst  case  24  hour  day. 
Given   the   following  variable   definitions, 

M(m,n)    -   the   number   of   days    in   the   nth  month   of    the  mth  year 

m  -  1   =  1980,    2   =  1981,    3    =  1982,    4  =  1983 

Y  -   the   number   of  years    in  the    total   period 

R  -   the  mean  ramping   rate   for   a   particular  hour 

D_.   -   summations   of    the    squared   24  hour   ramping   rates 

">,.   -  24  hour  mean  ramping   rates    (mean  24  hour   day) 

ayA   ~  24  hour   standard  deviations   of   ramping  rates 

S»  .   -   summations   of  24  hour   ramping  rates 

^24  -   total   number   of   days    in  period 

i  =  1.2,3 24 

the  means  and  standard  deviations  for  the  24  hour  rates  are 


. 


calculated   as   follows: 


1  12        M(L,k) 

111  hi 


iL 


L=l     k=l      j=l  S      (i) 

n,.(i)    =       =        MW/min  (3.4) 


24 


12  N24 


^         J       M(,l''n, 
m=l     n=l 


12        M(L,k) 

D      (i)    =        L         L         I  ftijkL  (MW/min)z  (3.5) 

L=l     k=l      j=l 


1         1         1 


D24(i)      -      (S24(i))2 

"24(i)    = 

I 

f                     N24  "  l 

MW/min  (3.6) 


In  addition,  the  two  worst  case  days  for  minimum  and  maximum 
ramping  rates  are  found  by  comparing  the  present  calculated  rate 
with  the  smallest  and  largest  rates  obseved  up  to  that  point.  If 
the  present  rate  is  smaller  or  larger  than  the  existing  m  in  or 
mar,  it  replaces  that  quantity  and  becomes  the  new  m  in  or  max. 
These  results  show  at  what  hours  of  the  day  worst  case  ramping 
occurred.  From  all  of  the  24  hour  results,  average  and  worst 
case  daily  ramping  trends  can  be  seen  and  compared  between  vari- 
ous power  production  schemes. 
12   Month  Ramping   Computations 

As    with    the    24   hour    ramping    computations,    the    12    month 
computations    create   a  year   in  which   each  month    is    represented   by 
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typical  and  worst  case  ramping  for  the  month.  Along  with  these 
quantities,  the  standard  deviation  of  all  ramping  for  a  particu- 
lar month-set  about  the  mean  monthly  ramp  is  also  calculated. 
The  term  month-set  emphasizes  the  fact  that  even  though  the  re- 
sults are  labeled  for  a  single  month,  up  to  three  months  of  data 
(three  months  with  the  same  name  for  three  years)  are  used  in 
calculations.  Also,  special  note  should  be  taken  to  the  mean  12 
month  results  which  incorporate  absolute  valued  ramping  rates. 
Since  the  average  ramping  rate  over  a  month  is  approximately  zero 
(as  many  positive  ramps  as  negative),  the  magnitude  of  ramping  is 
used,  thus  giving  information  on  absolute  average  monthly  trends. 
Given   the   following  variable   definitions, 

M(m,n)    -   the   number   of   days    in   the   nth  month   of    the  mth  year 

m  -  1   =  1980,    2   =  1981,    3   =  1982,    4  =  1983 

Y  -  the  number  of  years  in  the  total  period 

R  -  the  mean  ramping  rate  for  a  particular  hour 

V  -   summations  of  the  squared  12  month  ramping  rates 
m^,  -  12  month  mean  ramping  rates  (absolute) 

o\.2  ~  12  month  standard  deviations  of  ramping  rates 

S    -  summations  of  12  month  ramping  rates 

N-.  -  total  number  of  hours  in  month  period  (3  months) 

k  =  1,2.3,. ..,12 
the  means  and  standard  deviations  for  the  12  month  rates  are 
calculated  as  follows: 
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Y         M(L,k)    24 

1         1              1        'RijkLl 
L=l      j=l           i=l                          S      (k) 
■»,,<*>    =  =  MW/min        (3.7) 


"12  * 


N12(k) 


24    .    )         M(m,k) 
m=l 


Y         M(L,k)    24 

D12(k)    "       }         1  1  Ri.ikL  (MW/min,2 

L=l      j=l  i=l 


ijkL  (MW/min)  (3.8) 


D12(k)      -      (S12(k))2 


°12(k)      =      1/  N12(k)  MW/min  (3.9) 

N12(k)    -  1 

In  addition,  for  the  worst  case  minimum  and  maximum  12  month 
year,  a  similar  procedure  to  that  for  24  hour  calculations  is 
used.  Thes.e  results  yield  the  worst  ramping  which  occurred 
during  the  three-month  sets.  All  of  the  12  month  results  allow 
for  comparison  between  calculations  for  different  power  produc- 
tion schemes,  and  show  average  and  worst  case  monthly  trends. 
Weekend- Weekday  Ramp ing   Calculations 

To  conclude  the  ramping  computations,  weekend  and  weekday 
ramping  rates  are  calculated  to  give  results  that  compress  three 
years  of  ramping  rates  into  typical  and  worst  case  ramping  for  a 
weekend  and  a  weekday.  Also,  the  standard  deviations  for  the 
weekend-weekday  ramping  rates  are  calculated  with  their  respec- 
tive mean  value.  Again,  as  in  12  month  computations,  note  that 
the   mean  weekend-weekday   rates   are  formulated  using  absolute 
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ramping.      Given   the   following  variable   definitions, 

R  -   the  mean  ramping   rate   for   a   particular  honr 

D     -   summations   of    the    squared  weekend-weekday   ramping   rates 

m,   -  weekend-weekday  mean  ramping   rates    (absolute) 

c_   -  weekend-weekday  standard  deviations   of    ramping   rates 

S.    -   summations   of  weekend-weekday  ramping  rates 

N,   -   total    number   of   ramps    for  weekends    and  weekdays 

W  -  number   of  weeks    in  total   period 

s(k)   =     1     for     k=l  e(k)   =  2     for     k=l 

3  k=2  7  k=2 

k  =  1,2    (weekend,    weekday) 

the  means  and  standard  deviations  for   the  weekend-weekday  ramping 

rates   are   calculated  as   follows: 

W  e(k)         24 

L=l      j=s(k)    i=l  S    (k) 

"2(k)    =  "         MW/min      (3.10) 

24We(k)  N    (k) 


W  e(k)        24 

D2(U=  1         I  1  «?Jtt  (MW/min)2 

L-l     j=s(k)    i=l 


ijkL  (MW/min)^        (3.11) 


D2(k)      -      (S2(k))2 


°2(k)    =         II  N2(t)  MW/min  (3.12) 

N2(k)    -  1 

In  addition,    for    the  worst   case   weekend-weekday  ramping   computa- 
tions,   the    same   procedure    to   determine    these  values  as   used  above 
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is  implemented.  All  weekend-weekday  results  allow  for  compari- 
son between  calculations  for  different  power  production  schemes, 
and  show   average   and  worst  case  weekend-weekday  differences. 

Ramping  Distribution 

In  recaping  the  average  and  worst  case  ramping  calculations 
described  above,  it  was  noted  that  the  results  given  by  those 
calculations  allowed  perception  of  the  effect  that  added  wind 
generation  would  have  on  utility  operation.  However,  to  gain  an 
overall  picture  of  how  all  ramping  rates  within  a  given  range  are 
affected,  another  method  must  be  used.  This  method,  a  calculated 
ramping  distribution,  allows  insight  into  the  distributed  nature 
of  ramping  rates,  and  allows  closer  inspection  of  effects  due  to 
injected  wind  power   generation. 

Within  the  routines  for  this  thesis,  the  procedure  for 
finding  the  desired  distribution  involves  finding  the  probabilty 
density  function  (pdf)  of  the  ramping  rates,  and  then  calculating 
the  distribution  from  the  defining  equation  which  links  the  pdf 
and  the  probability  distribution.  As  seen  in  Figure  3.5,  an 
approximate  pdf  can  be  made  up  of  distinct  bins,  each  bin  being 
2s  wide  and  centered  about  o.  To  find  this  approximate  pdf,  a 
relation  mnst  be  found  between  the  probability  of  a  range  of 
ramps  occurring,  and  the  pdf  itself.  In  looking  at  Figure  3.5 
again,  note  that  for  proper  definition  of  probability  over  an 
interval   of   a   continuous   pdf,    that 
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..a+e 
P[a-e    <  R    <   a+e]    =  J  fR<r>dr  (3.13) 

a-e 

where  R  is  the  random  variable  for  ramping  rates.  However,  if 
the  bin  width  shown  in  Figure  3.5  is  small  enough,  then  the 
probability,  which  is  the  area  under  the  pdf  curve,  is  approxi- 
mately 


P[o-e    <  R   <  a+e]    =  2efR(a)  (3.14) 


Also,  for  large  N,  it  is  known  that 


_^   No.  of  r.  between  a-e  and  a+e 

P[a-e  <  R  <   a+e]  ■   * (3.15) 

N 


where  r  is  the  ith  discrete  ramp.  If  the  approximations  given 
by  equations  3.14  and  3.15  are  now  combined,  then  the  formula  for 
the  approximate  pdf  becomes 


^        No.  of  r.  between  a-e  and  a+e 

fR(o)  =  J (3.16) 

2N| 


where  fR(a)  is  the  approximate  pdf  of  the  ramping  rates  [6]. 

Now  that  the  pdf  has  been  defined  by  an  approximate  numeri- 
cal method,  the  distribution  can  be  found  using  the  relation 
between  the  pdf  and  the  distribution  function.  From  statistics, 
the  relation  between  continuous-probability  density  and  distribu- 
tion functions  is 
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FR(r)  =  f 


fR(a)da  (3.17) 
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where  fR(a)  is  the  value  of  the  pdf  at  a.  However,  this  distri- 
bution function  gives  the  probability  that  a  ramp  of  R  or  less 
has  been  observed  (integrating  from  0  to  r),  and  is  the  converse 
of  what  is  desired  (the  probability  of  observing  a  particular 
ramp  R      or   greater).      To  obtain  the  desired  result,    if  the  pdf 

were    integrated   from   r   to  »,    the   new   distribution    (F         (r))    would 

new 

be 


R(r)    =  1  -  j     fR(o)du     =  J     fR(a)da  (3"18) 


F        (r)   =  1  -  F 
new 


To  write  this  in  a  more  reasonable  form,  the  continuous  compo- 
nents must  be  replaced  with  their  respective  discrete  parts.  In 
substituting  2e  for  da,  fR(a)  (from  equation  3.16)  for  fR(a), 
summation  signs  for  integrals,  and  noting  that  the  2e's  cancell, 
equation   3.18  becomes 


_.     No.    of   r.    between   a-e   and   a+e 
Fnew(r)    =     I     * (3.19) 

N 


Now,  for  actual  calculations,  •  is  replaced  by  the  maximum  and 
minimum  values  of  +10.0  MW/min,  e  is  set  to  0.1  MW/min  for  a  0.2 
MW/min  bin  width,  r  is  restricted  to  the  range  -10.0  £  r  <.  10.0, 
and  a  is  incremented  by  +  0.2  MW/min.  Finally,  positive  and 
negative  ramping  distributions,  F  (r)  and  F„(r),  are  calculated 
separately,  and  are  represented  by 

10-0  v, 

^.    No.  of  r.  between  a-s  and  a+e 

Fp(r)  =  2 '0.0  <  r  <  10.0   (3.20) 

N 
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-10.0  „  . 

_         No.    of   r.    between  o-e   and   0+8 

FN(r)   =  1  ^~ ;-10.0  i  r  1  0.0    (3.21) 

(by  20?2) 


Earlier,  it  was  stated  that  within  routine  ealcnlations,  that  in 
calculating  the  distribution,  the  pdf  would  be  calculated  first. 
This  is  true  up  to  a  constant  value  divisor.  Instead  of  calcu- 
lating equation  3.16  as  it  is  written,  the  equation  times  2e  is 
calculated  (the  interior  of  equations  3.20  and  3.21)  and  thus 
allows  direct  summation  without  having  to  multiply  the  pdf  by  2e 
to    obtain    results. 

Upwards   Ramping  Control 

In  the  analysis  of  wind  power  affected  ramping  rates,  all 
analysis  so  far  has  assumed  that  turbine  operation  has  had  no 
outside  control.  In  actual  operation  though  this  will  most 
likely  not  be  the  case.  Within  this  section,  a  simple  control 
algorithm  that  limits  wind  turbine  upwards  ramping  by  spilling 
generated  power    is   discussed. 

As  shown  in  the  previous  section,  a  wind  turbine  produces 
power  with  respect  to  wind  speed  according  to  a  curve  like  that 
of  Figure  3.1,  with  a  resultant  wind  power  profile  as  seen  in 
Figure  3 J.  For  a  general  profile  of  mean  hourly  powers  (Figure 
3.6),  if  some  power  P2  is  greater  than  P  such  as  to  give  a 
ramping  rate  greater  than  some  desired  limit  RL,  then  the  control 
algorithm    spills  power   to  keep   that   ramping  rate  at  Rj  .      Spilling 
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of  wind  power  is  accomplished  mainly  through  blade  pitch  control. 
This  new  power  P2new  replaces  P  and  is  calcnlated  by 


2new 


Px  +  601^      MW  (3.22) 


where  R^  is  in  MW/min,  P  and  P  are  in  MW,  and  60  carries  the 
dimension  minutes.  Also,  to  measure  the  amount  of  lost  power  due 
to  this  ramping  control,  the  power  loss  is  summed  for  all  cases 
of  excess  ramping,  and  is  then  divided  by  the  total  amount  of 
power  production  for  the  period.  This  result  gives  the  percent 
power  loss  due  to  ramping  control,  and  allows  for  comparison 
between  different  power  production  schemes  and  different  amounts 
of   control. 

As  far  as  control  of  actual  wind  farms  is  concerned,  the 
technology  of  power  production  control  to  help  economic  dispatch, 
for  example,  has  been  in  use  for  at  least  20  years.  Thus,  the 
implementation  of  a  control  scheme  on  wind  power  production  would 
be  of  major  importance.  The  control  algorithm  used  in  this 
thesis  looks  only  at  ramping  rate  magnitudes,  and  does  not  take 
into  consideration  the  entire  power  system  as  a  whole.  Also,  for 
wind  arrays  consisting  of  two  or  more  farms,  the  control  algo- 
rithm used  in  this  thesis  monitors  the  aggregate  ramping  of  the 
wind  array  and  not  each  site  individually.  If  this  were  in 
actual  operation,  a  control  center  would  then  have  to  monitor 
each  individual  site,  add  their  individual  contributions  toge- 
ther, and  spill  power  at  the  most  convenient  site(s)  to  maintain 
economic   dispatch   and   keep   ramping   at  a   desired  level. 
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An  example  of  ramping  control  as  seen  on  a  wind  power  pro- 
file is  shown  in  Fignre  3.7.  The  control,  seen  as  positive 
parallel    slopes,     is    quite    evident    in   the   npper  plot. 

Correlation 

Besides  ramping  observations,  another  area  of  interest  is 
the  way  in  which  utility  load  demand  and  wind  power  production 
correlate.  That  is,  the  way  that  Kansas  load  demand  patterns 
after  Kansas  wind.  This  analysis  involves  autocorrelation  and 
crosscorrelation  functions,  and  treats  wind  and  utility  data  as 
discrete  signals  to  be  analyzed  in  a  signal  processing  method 
[6].  These  results  show  how  a  signal  correlates  with  itself,  and 
how    two    signals    correlate   together. 

For  the  analysis  performed  in  this  thesis,  the  procedure 
used  to  find  autocorrelations  is  shown  in  Figure  3.8,  while  the 
procedure  to  find  crosscorrelations  is  shown  in  Figure  3.9.  For 
both  procedures,  the  signals  are  preprocessed  by  subtracting  the 
mean  from  every  element  (see  Figure  3.10,  example  of  zero  mean 
utility  data),  or  taking  out  the  DC  component.  Next,  the  corre- 
lation is  performed  while  each  result  for  a  shift  k  is  divided  by 
the  product  of  the  standard  deviations  of  both  signals  (the 
product  for  the  autocorrelation  is  thus  the  square  of  the  single 
signal  standard  deviation).  This  has  the  effect  of  normalizing 
the  correlation  results  within  a  range  of  [+1.0,  -1.0].  Finally, 
note  that  the  correlation  functions  used  are  functions  of  time 
difference  only,  and  thus  show  how  correlation  varies  with  time 
shifting    of    data    sets. 
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Nutro    rui«z    3k-j-t->-    a«i--x     J 


In  equation  form,  the  left  (R    (k))  and  right   (i"    (k)) 

w»u  u,w 

side  crosscorrelations  are  written  as 

LeftSide  Rw  ^(k)   =     ^         w(i>u(i+k)  (3.23) 

V**  i=l 

Right  Side  Ru  w(k)    =     ]         u(i)w(i+k)  (3.24) 

Vw*  i-1 

where  u  is   the  otility  data  set,    w   is    the   wind  data   set,    N  is    the 
total  number  of   overlapping  valid  data  points,  «     is   the   standard 
deviation   of    the    utility   data,     and    a^    is    the    standard   deviation 
of   the  wind  data.      At   this  point,    the   most    confusing   argument    is 
that  of   the  difference  between  the   left-side  correlation  (denoted 
w,u)    and    the    right-side    correlation    (denoted   u,w).      If    the    util- 
ity  is   considered  the  base   by  which   the   equations  are  defined  for 
the    cross    correlation,    then   the    left    side    correlates    earlier 
occurring    winds    with   utility  data   of   a   later   time,    while    the 
right   side   correlates  later  occurring  winds  with  utility  data 
that   occurred  earlier.     For  example,    for   the  right  side   of   the 
crosscorrelation  plot    shown    in  Figure   3.11,    if    the   winds    for 
Tuttle    Creek   occurring    3    hours    later   blew    3   hours   earlier,     then 
the  utility  demand  and  wind  power  generation  would  correlate   on  a 
24  hour   cycle    (the   peak  close    to   the  center  of   the  plot   would 
indeed  be   on  zero   instead   of  being    shifted  by   three   hours). 

Finally,     to    calculate    the    autocorrelation   of    both    wind    and 
utility  data,    the   equation 
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Figure  3.11.  Utility  and  Tuttle  Creek  50m 
crosscorrelation. 
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Figure  3.12.   Plainville  50m  Autocorrelation. 
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1        <C 
Vx(k)   =     2     I      *<*>*<*+«  (3.25) 

A  1=1 

is  used,  where  i  is  the  data  being  analyzed,  a  is  the  standard 
deviation  of  the  data,  and  RJ>X  is  the  autocorrelation  of  the 
data.  In  showing  how  a  signal  correlates  with  itself,  important 
information  on  data  validity  can  be  extracted.  This  subject  will 
be  addressed  fully   in    the   next    section. 

An  example  plot  of  the  autocorrelation  for  Plainville  50  m 
is  shown  in  Figure  3.12.  As  can  be  seen  in  the  figure,  symmetry 
is  a  major  feature  of  the  autocorrlation  function.  Also,  in 
noting  that  the  horizontal  scale  is  in  hours,  this  plot  shows 
high  correlation  up   to  plus   or  minus    five  hours    time    shift. 

Data  Set  Validation 

Because  the  results  of  this  analysis  rely  heavily  on  the 
accuracy  of  the  wind  and  utility  data  used,  this  section  is 
included  to  verify  correctness  of  the  data  sets.  Both  the  utili- 
ty and  the  wind  data  suffer  from  problems  in  hourly  synchroniza- 
tion,   and   thus   both  data    sets   could  be   offset   up   to   two  hours. 

Since  the  start  of  wind  data  collection  at  Kansas  State 
University,  problems  with  data  transfer,  operator  error,  and 
transducer  failure  have  caused  an  incomplete  data  set  to  be 
collected.  Also,  when  data  was  available,  problems  in  synchroni- 
zation could  have  possibly  lead  to  offsets  of  up  to  two  hours 
between  consecntive  tapes  of  data.  Two  hours  is  chosen  as  the 
maximum   since   one  hour  definitely  arose   from   change   in    time. 


central  standard  to  central  daylight,  and  another  hour  is  added 
to  inclnde  any  other  possible  problems  that  might  have  occurred. 

Utility  data  on  the  other  hand  is  completely  synchronized 
with  its  respective  time  and  time  change.  However,  when  the 
utilities  make  the  time  change  from  standard  to  daylight,  move 
ahead  one  hour,  they  record  zero  energy  production  for  the  lost 
hour.  So,  instead  of  shifting  the  entire  data  set  back  to  cover 
this  hole,  an  interpolation  is  made  between  adjacent  hours  of 
energy  production,  and  the  situation  is  treated  as  though  there 
was  never  a  time  change.  At  first  this  would  not  seem  valid,  but 
since  the  utility  never  records  the  extra  energy  that  it  produces 
when  changing  from  daylight  to  standard  time,  the  hour  offset  due 
to  this  problem  only  occurs  for  at  most  six  months,  and  averages 
out  for  the  year.  Finally,  for  the  routines  used  in  this  analy- 
sis, it  is  important  to  note  that  the  utility  data  has  one  pro- 
jected hour  of  data  added  to  the  end  of  the  data  set.  This  is 
done  so  that  the  calculations  of  ramping  rates  can  extend  to  the 
last  available  hour  of  data. 

To  show  that  these  effects  will  not  alter  the  outcome  of 
results  significantly,  a  closer  look  at  the  autocorrelation  of 
wind  power  output  must  be  made.  In  looking  back  at  Figure  3.16 
it  is  found  that  wind  data  correlates  highly  (50  %  correlation) 
with  itself  even  when  shifted  in  time  by  up  to  5  hours,  and  for  a 
time  shift  of  2  hours,  is  even  better,  showing  80  %  correlation. 
Thus,  even  if  the  wind  data  were  offset  by  two  hours,  the  outcome 
would  not  be  significantly  different  since  the  major  variations 
of  wind  occur  beyond  a  time  shift  of  two  hours.   In  fact,  from 
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Figure  3.16,  major  variations  (periodic  trends)  appear  to  be 
occurring  on  24  hour  and  4  day  cycles.  So,  since  the  utility 
offset  averages  out  over  a  yearly  period  and  wind  correlates  well 
with  itself  up  to  a  2  hour  time  shift,  results  obtained  froa 
these    two  data   sets   will   express   true   occurrences. 
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CHAPTER  IV 
RESULTS  AND  DISCUSSION 

After  mating  numerous  runs  of  each  of  the  main  routines, 
results  for  ramping,  and  correlation  were  compiled  and  analyzed. 
These  results  constitute  the  heart  of  this  thesis,  and  show  the 
effects  that  added  wind  generation  could  have  on  utility  opera- 
tion. 

For  all  of  the  results  discussed  below,  standard  input 
parameters  for  the  wind  power  calculations  will  be  used  unless 
otherwise  noted.  These  parameters  consist  of  a  300  MW  rated  wind 
array  power  output,  a  cut-in  turbine  wind  speed  of  6.0  mps,  a 
rated  turbine  wind  speed  of  12.0  mps,  no  upwards  ramping  control, 
and  a  wind  array  of  one   farm. 

Utility  Ramping 

Utility  ramping  by  itself  is  needed  as  a  base  by  which 
comparisons  can  be  made.  In  following  analyses,  this  base  result 
is  often  termed  the  utility  norm  or  unaffected  utility,  and 
represents  the  normal  case  of  utility  operation  without  any  type 
of  added  wind  generation. 
24  Hour  Results 

With  the  utility  24  hour  ramping  rates,  the  first  phenomena 
seen  should  be  that  of  the  ever  present  morning  pickup  of  Kansas 
load  demand.  In  looking  at  Figure  4.1  and  noting  that  the  hori- 
zontal axis  increments  are  1.15  hours,  it  is  found  that  the 
morning  pickup,    with  its      highly  visible  first  positive  hump, 
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appears  to  have  an  average  peak  of  2.63  MW/min.  Also  seen  in 
this  plot  are  the  noon  and  evening  pickups  (second  and  third 
positive  humps),  while  the  evening  fall  off  obtains  an  average 
negative  peak  of  -2.88  MW/min.  Another  interesting  result  can  be 
observed  in  the  plot  of  the  standard  deviation.  As  can  be  seen 
in  Figure  4.1,  the  greatest  variation  in  ramping  occurs  during 
the  morning  pickup  between  6  and  7  am  (some  are  contemplating 
getting  out  of  bed  for  work),  while  the  least  amount  of  variation 
occurs  during  the  tail  end  of  the  work  day  (everyone  takes  a 
break  or  goes  home  from  work)  and  during  the  early  morning  hours 
(everyone  is  asleep).  These  results  are  not  in  the  least  sur- 
prising, but  are  none  the  less  humorous.  For  maximum  ramping,  it 
is  also  seen  that  worst  case  situations  occur  during  the  morning 
pickup  with  a  peak  of  5.43  MW/min.  Since  this  peak  occurs 
around  8  am,  a  suspective  cause  would  be  that  of  the  start  to  the 
eight  hour  work  day.  Finally,  the  minimum  ramping  calculations 
show  worst  case  fall  off  to  occur  at  10  pm  and  have  a  minimum 
value  of  -5.95  MW/min.  A  suspective  cause  again  is  that  of 
industrial  nature. 
12  Month  Ramping   Results 

In  looking  at  ramping  rates  on  a  monthly  time  axis,  intui- 
tion suggests  that  the  number  of  variations  will  decrease.  This 
would  be  due  mainly  to  the  averaging  out  of  the  diurnal  and 
weekly  cycles,  thus  showing  only  seasonal  variations.  The  12 
month  results  in  Figure  4.2  have  this  characteristic,  and  show 
one  major  trend,    that  being  all  around  increases  during  summer 
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months.  Thus,  during  this  time  there  is  increased  ramping  and 
increased  deviation,  while  there  are  lower  minimums  and  higher 
maximums.  The  worst  case  values  for  the  minimums  and  maximums 
match  those  of  the  worst  case  24  hour  results,  and  occur  during 
the  evening  fall  off  and  morning  pickup  respectively. 
Weekend-Weekday  Results 

Concerning  these  results,  there  appears  to  be  hardly  any 
noticeable  differences  between  the  weekdays  and  the  weekends. 
However,  as  seen  in  Figure  4.3,  the  weekends  do  show  a  slight  in- 
crease in  both  the  mean  value  and  the  deviation  from  that  mean 
value,  while  the  weekdays  show  slightly  higher  worst  case  ramping 
conditions.  This  matches  what  would  be  expected,  since  during 
the  weekend  the  8  hour  work  day  pattern  is  not  as  prevalent 
causing  more  random  and  spontaneous  power  demand.  Also,  the 
greater  minimum  and  maximum  ramping  occurrence  during  the  week 
would  seem  to  stem  from  industrial  operations  which  occur  during 
this  time. 
Ramping  Distribution  Results 

As  in  the  12  month  results  just  mentioned  above,  the  distri- 
bution (Figure  4.4)  has  only  one  major  trend,  that  being  exponen- 
tial in  nature,  and  decreasing  towards  higher  ramping  rates. 
This  trend  shows  that  higher  ramping  rates  have  a  decreasing 
exponential  probability  of  occurrence  and  indicates  utility  ramp- 
ing rarely   to  exceed   5.0    MW/min. 

Wind  Power  Ramping 

Wind  power  ramping  results  are  used  primarily  to  show  the 
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WEEKEND-WEEKDAY  MEANS  (ABSOLUTE) 
WEEKEND  ->    1.22  WEEKDAY  - 


1.11 


WEEKEND-WEEKDAY  STANDARD  DEVIATIONS  (ABSOLUTE) 
WEEKEND  — >    1,0?  WEEKDAY  — >     .9 

WEEKEND-WEEKDAY  MINIMUMS 
WEEKEND  — >   -5.93  WEEKDAY  -->   -5.? 

WEEKEND-WEEKDAY  MAXIMUMS 
WEEKEND  — >    5.17  WEEKDAY  — >    5,4 


Figure  4.3.  Weekend-weekday  results  for  utility 
power  production. 
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ramping  effects  of  wind  power  production  alone.  Since  wind  is 
fairly  random,  one  would  tend  to  think  the  ramping  results  from 
such  power  production  to  be  random  also.  However,  to  take  a 
complete  look  at  the  wind  and  verify  such  ideas,  the  following 
analysis  of  wind  power  ramping  is  made. 
24  Hour  Results 

By  inspecting  the  24  hour  results  for  all  the  sites,  it  was 
found  that  wind  ramping  rates,  as  wind  itself,  are  very  random. 
As  seen  in  Figure  4.5  for  Atlanta  50  m,  the  mean  ramping  rates 
have  no  apparent  pattern,  while  the  standard  deviations,  mini- 
mums,  and  maximums  do  the  same.  From  the  minimums  and  maximums 
however,  the  worst  case  ramping  can  be  averaged,  and  is  observed 
to  be  around  -3.0  MW/min  and  3.0  MW/min  respectively.  Even 
though  only  one  site  is  shown  in  the  figure  for  24  hour  ramping 
computations,  all  other  sites  have  approximately  the  same  re- 
sults. 
12  Month  Ramping   Results 

For  the  12  month  ramping  rates,  the  diurnal  cycle  is  aver- 
aged out  again,  thus  allowing  seasonal  variations  to  show.  As 
seen  in  Figure  4.6  for  Tuttle  Creek  50  m  however,  the  mean  12 
month  rates  show  no  seasonal  variation  whatsoever,  and  appear  to 
be  rather  uniform  in  nature  over  the  entire  yearly  period.  The 
standard  deviation  also  has  this  character,  but  shows  a  slight 
increase  in  deviation  towards  the  center  of  the  year.  Finally  in 
the  minimums  and  maximums,  some  trends  begin  to  show.  As  seen  in 
the  minumum   ramping  plot,    the  early  summer  months   show   the  worst 
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sustained  negative  ramping  conditions  for  the  entire  yeax,  reach- 
ing -5.0  MW/min  (due  to  drop  from  rated  power  of  300  MW  to  zero 
power  in  one  hour).  Also,  the  postitive  rates  show  this  effect 
for  the  early  summer  months,  and  reach  the  opposite  extreme  of 
5.0  MW/min  (dne  to  increase  from  zero  power  to  rated  power  of  300 
MW  in  one  hour).  So,  even  though  wind  power  is  a  randomly  occur- 
ring variable,  there  is  some  seasonal  variation  in  its  produc- 
tion. 

Utility  Itamping  After  Wind  Power  Injection 

To  analyze  resultant  utility  ramping  effects,  three  of  the 
major  calculations  used  above  and  control  of  wind  power  produc- 
tion parameters  will  be  nsed.  The  calculations  consist  of  24 
hour  ramping  rates,  ramping  distributions,  and  percent  energy 
loss,  while  control  of  the  wind  power  production  will  be  accom- 
plished through  siting  of  a  wind  farm,  choosing  cut-in  and  rated 
wind  speeds  for  a  turbine,  allowing  specific  amounts  of  wind 
array  penetration,  and  increasing  the  number  of  farms  in  a  wind 
"ray-  Note  for  the  discussion  below,  that  12  month  and  weekend- 
weekday  results  will  not  be  used.  This  is  done  since  12  month 
results  are  not  complete  for  all  of  the  sites,  due  to  incomplete 
data  sets,  and  weekend-weekday  results  for  added  wind  power 
show  random  effects. 
Siting   Control 

Siting  control  involves  taking  a  look  at  the  results  from 
different  sites  and  determining  the  optimum  site  for  location  of 
a   farm,    within  the   site    limitations  of   this   study.      At  first,    one 
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would  think  of  this  as  applying  no  control  to  the  injected  wind 
power,  but  when  looked  at  on  a  broader  basis  of  site  to  site,  is 
actually  control  of  a  specification  for  implemented  wind  genera- 
tion. 

Referring  to  Figure  4.7,  the  ramping  distribution  for  the  50 
m  results  tend  to  be  exponentially  decreasing  towards  higher 
ramping  conditions.  Tuttle  Creek  and  Plainville  create  increased 
utility  ramping  over  the  entire  range  of  rates,  and  for  2.0 
MW/min  or  greater,  both  show  approximately  a  10%  increase  in  the 
probability  of  occurrence  for  negative  rates.  Atlanta  affects 
the  utility  only  lightly  in  the  positive  rates,  and  for  2.0 
MW/min  or  greater,  begins  to  help  the  utility  with  negative 
ramping.  This  feature  of  added  wind  generation  aiding  utility 
operation  is  either  a  statistical  mishap,  in  that  all  of  the  data 
for  this  particular  time  period  gave  such  a  result,  or  Atlanta  is 
truly  a  beneficial    site   for    the   production   of   wind   energy. 

For  the  30  m  distribution  results,  all  sites  raise  the 
utility's  distributed  ramping  rates  for  both  positive  and  nega- 
tive ramping.  Although  the  curves  are  jumbled  together,  Tuttle 
Creek  shows  an  increased  negative  effect  on  the  utility  beyond 
that  of  the  other  three  sites  for  negative  ramping.  However, 
since  this  difference  is  very  small,  the  four  sites  at  30  m 
show   signs   of  having   coincident   effects   on  utility  ramping. 

In  continuing  observance  of  site  effects,  if  the  24  hour 
mean  ramping  rates  are  inspected  (Figure  4.8),  of  the  three  sites 
at    50    m,    Atlanta    again    shows    to    be    the    most    varied    from    the 


46 


47 


r 

o 

c 

u. 

r 

£ 

a 

o 

S 

£ 

z 

1- 

to 

S 
IB 

» 

Y 

E 

Lll 

Ui 

;» 

111 

III 

2 

HI 

£C 

-1 

a 

a 

U 

J 

N 

« 

» 

n 

UI 

> 

1- 

2 

J 

Z 

Z 

"! 


vS 


(0 

HI 

I- 

s 
e 


z 


ojifcoiiiZH  oil  i-xai  i-«aij 


Q-UjauiUZl-    OIL    KTLli    (-KEU1 


Figore  4.7.  Dtility  ramping  distribution  after  wind  power 
injection. 
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Figure  4.8.  Utility  24  hour  means  after  wind 


power  injection. 


utility  norm.  Atlanta  tends  to  make  positive  ramping  increasingly 
positive,  and  varies  on  its  effects  to  negative  ramping.  Plain- 
ville  and  Tuttle  Creek  also  canse  tie  utility  to  stray  away  from 
average  operation,  but  have  decreased  effects  as  compared  to 
those  of  Atlanta.  In  fact,  Tuttle  Creek  at  50  m  has  very  little 
impact  on  the  average  operation  of  the  utility  as  seen  in  the 
plot  of  24  hour  mean  results. 

The  30  m  24  hour  mean  results  compare  directly  to  the  dis- 
tributions for  30  m  towers  in  that  they  all  appear  to  act  the 
same.  This  time,  Wright  pulls  away  from  the  group  by  having  a 
slightly  higher  peak  ramp  during  the  evening  dinner  hour,  and  by 
altering  the  shape  of  the  morning  pickup.  Still  though,  the  30  n 
effects  show  to  be  most  consistent  between  sites. 

Moving  to  the  standard  deviations  of  the  24  hour  results 
(Figure  4.9),  in  the  50  m  case,  Atlanta  again  separates  from 
Tuttle  Creek  and  Plainville.  Instead  of  having  a  complete  in- 
crease in  standard  deviation,  Atlanta's  effect  on  the  utility 
causes  the  standard  deviation  to  fall  below  that  of  the  utility 
norm  during  the  hours  of  7  am  to  1  pa.  Tuttle  Creek  and  Plain- 
ville however,  cause  ramping  rates  to  deviate  greater  during  this 
time  by  0.2  MW/min,  and  cause  the  utility  to  maintain  higher 
deviations  for  the  entire  24  hour  period.  This  is  also  the  case 
for  all  30  m  standard  deviations,  with  the  exception  of  Plain- 
ville which  dips  below  the  unaffected  utility  for  one  hour. 

For  the  worst  case  minimums  and  maximums,  all  sites,  both  30 
and  50  m  (Figures  4.10  and  4.11),  cause  greater  extreme  ramping 
cases  for  most  of  their  24  hour  periods.   Also,  an  interesting 
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Figure  4.9.  Utility  24  hour  standard  deviations  after  wind 
power  injection. 
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Figure  4.10.  Utility  24  hour  minimums  after  wind  power  injection. 
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Figure  4.11.  Utility  24  hour  maximums  after  wind  power  injection. 
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outcome  is  seen  in  Atlanta's  effect  on  the  utility  for  the  24 
hour  maximums.  Atlanta  causes  utility  maximums  to  decrease  again 
from  7  am  to  1  pi,  but  during  the  morning  and  evening  pickups, 
causes  a  higher  extreme  ramping  case  to  occur,  above  even  those 
of  Plainville  and  Tuttle  Creek.  So,  as  seen  earlier,  Atlanta  was 
an  excellent  site  when  the  distributions  were  studied,  but  shows 
to  have  the  greatest  extreme  ramping.  This  result  is  very  possi- 
ble since  the  distribution  would  most  likely  not  register  the  few 
worst  case  conditions  occurring  at  higher  ramping  rates. 
Cut-in  and  Rated  Wind  Speed   Control 

For  the  cut-in  and  rated  wind  speed  control  results  and  all 
results  to  follow,  a  base  core  of  sites  and  tower  heights  will  be 
used.  This  is  done  to  limit  the  number  of  test  cases  that  must 
be  run  to  analyze  the  results  of  wind  power  injection  control. 
The  sites  to  be  used  will  be  Tuttle  Creek  50  m,  Plainville  50  m, 
Atlanta  50  m,  and  Wright  30  m.  It  is  important  to  note  that  only 
one  30  m  tower  will  be  used,  since  all  30  m  towers  for  the 
various  locations  have  relatively  the  same  ramping  characteris- 
tics   as    discused    in   the   previous    section. 

Starting  with  the  ramping  distributions,  it  is  seen  that 
increasing  uc  and  ur  causes  a  proportional  decrease  in  the  ramp- 
ing distributions  for  all  sites  (Figure  4.12).  It  would  seem 
from  this  that  higher  u0  and  ur  perform  a  filtering  action  for 
low  wind  speed  cut-ins  (which  create  a  higher  noise  component). 
Atlanta  (Figure  4.13)  shows  the  same  trend  of  a  lower  distribu- 
tion for  higher  uo   and  ur,    but   as    seen  earlier,    drops    the   utility 
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Figure   4.12.    Utility    ramping   distributions   after  Tuttle   Creek 
injection  for    three  valnes    of   u     and  u   . 
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Figure   4.13.    Utility  ramping   distribution  after  Atlanta    injection 
for   three  values   of  u     and  u    . 
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distribution  below    that   of   the   unaffected  ramping   case. 

For  the  average  24  hour  day,  changes  in  u  and  n  do  not 
create  any  visible  trend.  Injected  wind  power  does  alter  the 
shape  of  the  mean  24  hour  curve,  but  as  u  and  u  are  varied,  the 
curves  show  no  deviation  from  the  standard  parameter  curve  dis- 
cussed under  site  control.  However,  the  standard  deviation  to 
this  mean  curve  (Figures  4.14  and  4.15)  shows  a  definite  trend  of 
increased  ramping  deviation  due  to  decreases  in  u  and  u  .  This 
trend  is  seen  for  all  sites,  and  in  the  worst  case  for  Tuttle 
Creek,  increases  the  standard  deviation  by  0.5  MW/min  during  the 
morning    pickup. 

For  the  24  hour  minimums  and  maximums,  all  sites  show  the 
majority  of  ramping  to  be  more  severe.  No  real  trend  can  be 
spotted  for  variations  in  uc  and  ur,  while  the  injected  wind 
power  again  shows  its  negative  affect  upon  the  utility  worst  ease 
ramping  conditions.  An  interesting  effect  is  also  noted  at 
Atlanta,  where  the  greatest  worst  case  ramp  of  8.5  MW/min  oc- 
cured.  All  other  sites  stayed  well  below  the  8.0  MW/min  level. 
Penetration  Control 

As  with  the  variation  of  uc  and  ur,  the  variation  of  pene- 
tration causes  a  dramatic  trend  to  occur  on  the  distribution  of 
ramping  rates.  For  all  sites,  as  wind  power  penetration  is 
increased,  both  negative  and  positive  ramping  distributions  rise 
as  a  result.  In  fact,  for  the  Tuttle  Creek  (Figures  4.16  and 
4.17),  Plainville,  and  Wright  positive  distributions,  the  change 
from  zero  to  15%  wind  penetration  causes  a  10%  increase  in  obser- 
vation of   ramping  rates    greater   than  or    equal    to   2.0    MW/min,    and 
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Figure  4.14.  Utility  24  hour  standard  deviations  after 
Creek  injection  for  three  values  of  a  and 
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Figure  4.15.    Utility  24  hour   standard  deviations   after  Atlanta 
injection  for    three  values   of   u     and  u   . 
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Figure  4.16.  Utility  ramping  distributions  after  Tattle  Creek 
injection  at  three  penetration  levels. 
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Figure   4.17.    Blow  up   of  Figure  4.16. 


likewise  for  the  negative  distribution.  Atlanta  again  has  the 
feature  of  falling  below  the  unaffected  utility  distribution,  but 
for  the  positive  distribution,  moves  entirely  above  for  3.0 
MW/min  or  greater. 

For  the  mean  24  hour  day,  all  sites  show  no  trend  for  in- 
creasing penetration.  This  would  tend  to  show  that  the  average 
value  of  wind  power  as  calculated  by  24  hour  periodic  sampling, 
is  a  constant,  and  is  not  dependent  on  the  amplitude  of  the  power 
production.  From  the  24  hour  plots,  note  that  even  though  there 
are  no  noticeable  variations  with  respect  to  penetration,  the 
injected  wind  power  still  causes  the  utility  to  deviate  from  its 
standard  24  hour   mean   profile. 

As  expected,  increasing  penetration,  like  the  distributions, 
causes  all  around  increases  in  the  24  hour  standard  deviation, 
with  Atlanta  and  a  few  other  rare  cases  showing  values  below  the 
utility  norm.  As  seen  at  Tuttle  Creek,  an  increase  from  zero  to 
15%  penetration  causes  the  standard  deviation  to  increase  by  0.5 
MW/min.  Also,  the  minimums  and  maximums  (Figure  4.19)  had  the 
same  increasing  trend,  but  showed  most  notice  at  the  higher 
penetration  levels  of  10  and  15%.  As  an  example,  during  the 
morning  pickup  for  Atlanta,  worst  case  ramping  increased  from  the 
utility  norm  of  5.5  MW/min  to  6.0  MW/min  at  5%  penetration,  8.25 
MW/min  at  10%  penetration,  and  finally  to  just  under  11.0  MW/min 
at  15%  penetration.  For  minimum  ramping,  Tuttle  Creek  (Figure 
4.18)  had  the  worst  impact  causing  a  -10.5  MW/min  ramp  during  the 
late    evening    fall    off. 


61 


62 


Figure  4.18.  Utility  24  hour  minimums  after  Tuttle  Cre 
injection  at  three  penetration  levels 
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Figure  4.19.  Utility  24  tour  maximums  after  Tuttle  Creek 
injection  at  three  penetration  levels. 


Upwards   Ramp inn   Control 

As  noted  in  the  procedures  section,  resultant  utility  power 
production  is  calculated  by  subtracting  wind  power  from  the 
Kansas  load  demand.  Thus ,  when  upwards  wind  power  ramping  is 
controlled,  the  main  effect  appears  on  cases  of  negative  utility 
ramping.  That  is,  the  utility  will  decrease  power  output  to  allow 
wind  power  to  assume  more  of  the  load.  So,  in  looking  at  the 
distributions  of  the  modified  utility,  inspection  shows  that 
indeed  the  postive  ramping  distribution  remains  unchanged,  while 
the  negative  distribution  shows  a  decreased  trend  with  increasing 
amounts  of  control.  The  numerical  value  of  control  in  MW/min 
will  actually  get  less,  but  represents  increased  control.  This 
trend  (Figure  4.20)  is  observed  for  all  sites,  while  Atlanta 
again  begins  to  dip  below  the  negative  distribution  for  the 
utility  norm  at   2.0   MW/min. 

As  with  the  penetration  variations,  changes  in  control  do 
not  create  trends  in  the  average  24  hour  day.  These  successive 
results  would  tend  to  show  that  the  mean  24  hour  day  is  a  charac- 
teristic mainly  of  the  site  and  not  of  wind  power  production 
methods.  However,  as  seen  in  the  standard  deviation  for  the  24 
hour  results,  increased  control  helps  bring  widely  deviating  wind 
induced  ramps  closer  to  those  of  the  utility  norm.  This  trend  is 
seen  for  all  sites,  while  Atlanta  still  shows  signs  of  aiding 
utility  ramping   deviations   from    morning  till  noon. 

For  the  minimum  24  hour  ramps,  a  very  beneficial  trend 
develops.  As  control  is  increased,  the  worst  case  minimums  of 
the  utility  affected  by  wind  power   begin   to   approach    those   of    the 
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Figure  4.20.  Negative  utility  ramping  distribution  after  Tuttle 
Creek  injection  at  three  penetration  levels. 


utility  norm  (Figure  4.21).  In  fact,  for  tie  0.5  MW/min  control, 
all  sites  aid  the  utility  in  worst  case  ramps  for  particular 
hours.  This  would  tend  to  be  so  since  wind  power  that  picks  up 
and  adds  to  meet  Kansas  load  demand,  should  cause  the  utility  to 
ramp  with  less  intensity  as  long  as  the  wind  power  is  added  at  a 
steady  rate.  Another  important  effect  that  control  has  on  the 
minimums  is  that  of  decreasing  the  spiked  behavior  of  the  worst 
case  ramps  of  the  affected  utility.  In  essence,  the  control  acts 
as  a  low  ramping  rate  pass  filter  for  negative  ramping.  The  24 
hour  maximums  however,  show  no  evident  trends  for  change  in 
control,  and  only  in  a  few  instances  do  the  sites  create  slight 
decreases    in  worst   case  ramping. 

Now  that  the  ramping  characteristics  for  ramping  control 
have  been  addressed,  a  look  at  the  percent  energy  loss  due  to 
upwards  ramping  control  can  be  made.  For  each  site  at  50  m 
(Figure  4.22),  it  can  be  seen  that  this  percent  loss  decreases 
from  Tuttle  Creek,  to  Plainville,  to  Atlanta,  while  at  30  m 
decreases  from  Tuttle  Creek,  to  Plainville,  to  Wright,  to  Atlan- 
ta. For  0.5  MW/min  control.  Tuttle  Creek  shows  the  greatest 
percent  energy  loss  of  18.0%  at  50  m.  while  Wright  30  m  gives  the 
least  amount  of  loss  at  11.4%. 
ffiud  Array  Control    (multiple   farm) 

As  seen  with  upwards  ramping  control,  the  effect  of  added 
wind  generation  can  be  smoothed  by  limiting  wind  turbine  ramping. 
Another  method  used  to  solve  this  problem,  is  that  of  multiple 
farm   wind   arrays,    and   for    this    discussion,    will   be    accomplished 
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Figure  4.21.  Utility  24  hour  minimums  after  Tuttle  Creek 
injection  with  three  levels  of  control. 
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Figure   4.22.   Percent    energy   loss   dne    to   control    for   50m  sites. 
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by  using  up  to  three  sites  for  wind  power  production.  The  resul- 
tant effects  caused  by  this  injected  wind  power  will  be  calcu- 
lated and  analyzed  for  a  single  site  wind  array  at  Tutfle  Creek 
30  n>,  a  two  site  wind  array  consisting  of  Tuttle  Creek  and  Wright 
30  m,  and  a  three  site  wind  array  consisting  of  Tuttle  Creek, 
Wright,  and  Plainville  30  m.  Note  that  30  m  sites  are  used  in 
this  analysis  since  no  three  site  overlap  is  available  for  valid 
SO  m  data. 

Proceeding  to  the  distribution  results  (Figures  4.23  and 
4.24),  the  two  site  negative  distribution  begins  to  align  with 
the  utility  norm  around  3.75  MW/min,  while  the  three  site  wind 
array  actually  falls  below  at  this  point,  producing  a  much  de- 
sired result.  For  the  positive  rates,  at  first  the  two  and  three 
site  wind  array  effects  follow  the  distribution  of  the  single 
site  effect,  but  fall  off  rapidly  around  4.5  MW/min,  and  finally 
approach  the  unaffected  utility  distribution.  Also,  the  three 
site  wind  array  appears  to  fall  to  the  utility  norm  much  quicker 
than  the  two  site  array,  thus  showing  that  increasing  the  number 
of    sites    decreases   the    likelihood  of   experiencing  higher   ramps. 

For  the  24  hour  means,  as  more  sites  are  combined,  aggregate 
Kansas  wind  regime  characteristics  become  visible.  These  varia- 
tions consist  of  minor  modifications  to  the  morning  pickup,  with 
a  major  increase  in  utility  ramping  needed  during  the  evening 
pickup,    when  overall   winds   tend   to   decrease. 

As  seen  earlier  for  ramping  control,  standard  deviations  of 
the  24  hour  rates  show  significant  decreases  for  two  and  three 
site    wind    arrays.       For    about    five    hours    out    of    the    day,     the 
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Figure   4.23.   Utility  ramping   distributions    after   injected  wind 
power   from  1,    2,    and    3   site  wind  arrays. 
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Figure  4.24.    Blow  up  of   Figure   4.23 
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deviation  is  shown  to  fall  just  below  that  of  the  unaffected 
utility,  and  in  only  one  case,  the  evening  pickup,  increases  the 
standard  deviation.  This  increase  in  deviation  would  tend  to 
support  the  idea  that  in  terms  of  the  aggregate  Kansas  wind 
regime,  winds  fall  off  fairly  rapidly  during  the  early  evening 
hours. 

Proceeding  to  the  24  hour  minimums  (Figure  4.25),  almost  all 
hours  show  aid  to  the  utility  for  three  site  injection,  while 
increasing  the  number  of  sites  greatly  decreases  worst  case 
effects  due  to  a  single  site  injection.  This  also  happens  with 
the  maximums  (Figure  4.26),  but  aid  to  the  utility  only  occurs 
about  half  the  time  for  the  three  site  array,  and  one  third  of 
the  time  for  the  two  site  array.  The  two  and  three  site  wind 
arrays  also  smooth  the  spiked  behavior  of  worst  case  ramping 
conditions    as    seen   for    a   single    site    injection. 

In  terms  of  percent  energy  loss,  as  the  number  of  sites  in 
an  array  increases,  the  energy  loss  due  to  upwards  ramping  con- 
trol decreases.  As  seen  on  the  plot  of  percent  energy  lost 
versus  control  (Figure  4.27)  for  a  three  site  array,  percent  loss 
drops  below  5.0%  for  0.5  MW/min  control,  and  below  1.0%  for  1 
MW/min  control.  So,  as  wind  power  ramping  becomes  smoothed  by 
multiple  farms,  less  power  is  lost  to  control. 
Aggregate   Control   Example 

This  example  shows  that  if  beneficial  parameters  for  sites 
(wind  array  size,  uo  and  ur,  ramping  control)  are  chosen,  that 
higher   penetration   effects    can  be  made    to  match   those    of    a    lower 


Figure  4.25.  Utility  24  hour  minimums  after  injected  wind 
from  1,  2,  and  3  site  wind  arrays. 
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Figure  4.26.  Utility  24  hour  maximums  after  injected  wind  power 
from  1,  2,  and  3  site  wind  arrays. 
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Figure   4.27.   Percent   energy    loss   due    to   control   on  1.    2,    and   3 
site  wind  arrays. 
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penetration  level.  To  show  this  effect,  the  control  model  will 
consist  of  a  three  site  wind  array  containing  Tnttle  Creek, 
Wright,  and  Plainville  30  m,  all  with  uc  =  7  mps  and  n  =  14  mps, 
at  15%  wind  array  penetration,  and  finally  having  1.0  MW/min 
ramping  control.  The  results  from  this  case  will  be  compared 
against    three    other    test    cases: 

1.  Tnttle   Creek  30  m,    uc   =   5  mps,    u     =  10  mps,    no   control, 
and    5%  penetration. 

2.  Utility  without   any   added  wind   generation. 

3.  Tnttle   Creek  30  m,    u0    ■  7  mps,    u     =  14  mps,    no   control, 
and   15%  penetration. 

Beginning  with  the  positive  ramping  distribution,  it  is  seen 
that  the  three  site  smoothing  effect  starts  to  occur  around  4.0 
to  5.0  MW/min,  and  causes  the  affected  distribution  to  approach 
the  utility  norm  and  the  5%  penetration  case  before  reaching  5.2 
MW/min.  Earlier  in  the  distribution,  the  control  model  was  well 
above  all  other  cases,  but  proves  to  be  satisfactory  as  long  as 
it    decreases    before   reaching  higher   ramping   conditions. 

For  the  negative  distribution,  the  control  example  homes  in 
much  quicker  to  the  utility  norm  around  4.0  MW/min,  and  shows  to 
be  a  beneficial  power  source  after  this  point.  This  trend  was 
also  seen  under  the  ramping  control,  and  so  must  be  a  combined 
smoothing  effect  from  the  three  site  array  and  upwards  ramping 
control. 

As  seen  earlier,  the  24  hour  means  show  the  difference 
between  a  single  site  characteristic  and  that  of  the  aggregate 
Kansas  wind  regime.  Also,  aggregate  injection  again  causes  the 
utility   to    deviate    further    from    the    utility   norm    during    evening 
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pickup.  Likewise,  the  standard  deviation  for  the  control  model 
shows  the  aggregate  effect  by  going  above  all  cases  during  the 
evening  pickup.  However,  for  most  hours  the  control  model  shows 
aid  to  the  utility,  and  results  in  following  the  5%  penetration 
case    closely. 

Continuing  with  24  hour  minimums  and  maximums,  it  is  found 
that  both  are  brought  under  control.  In  fact,  the  control  model 
shows  aid  to  the  utility  minimums  for  over  half  of  the  worst  case 
24  hour  day.  Thus,  the  minimums  have  now  been  brought  down  to  a 
reasonable  level  which  can  be  dealt  with  by  the  utility.  The 
maximums  (Figure  4.28)  are  affected  in  the  same  way  as  the  mini- 
mums, showing  only  two  cases  of  going  slightly  above  the  5% 
penetration  maximums. 

Finally,  an  important  note  shoud  be  made  regarding  upwards 
ramping  control.  Since  this  type  of  control  was  incorporated, 
there  will  be  a  detremental  effect  to  the  efficiency  of  produced 
wind  generation.  However,  when  all  power  loss  due  to  upwards 
ramping  control  is  summed  and  divided  by  the  total  power  produc- 
tion, it  is  found  that  for  this  example  with  1.0  MW/min  ramping 
control,    only  2.01%  of    the    total    energy  production  was   lost. 

Correlation 

To  conclude  the  Results  and  Discussion,  a  study  of  wind  and 
utility  power  correlation  will  be  addressed.  First,  utility  and 
wind  power  autocorrelations  will  be  discussed  to  observe  the 
components  of  correlation  that  each  has  individually.   Then,  the 
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Figure  4.28.  Utility  24  hour  maximums  after  injection  of  wind 
power  from  aggregate  control  example. 
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crosscorrelation  of  the  utility  power  with  various  sites  of 
standard  power  production  parameters  will  be  analyzed  for  trends 
and  patterns.  All  of  the  correlation  plots  generated  for  this 
analysis  have  a. maximum  time  shift  of  one  week  (168  hours)  in 
both  the  positive  and  the  negative  direction.  This  will  not 
allow  detection  of  any  trends  that  might  occur  beyond  a  one  week 
time  shift. 
Autocorrelation 

Beginning  with  the  utility  autocorrelation,  it  is  seen  that 
all  correlation  values  for  increasing  time  shift  are  positive. 
This  is  an  interesting  phenomena,  but  turns  out  to  be  a  trivial 
subsequence  of  the  zero  mean  data  set.  Since  the  time  shift  for 
the  autocorrelation  has  a  maximum  shift  of  one  week,  the  values 
will  stay  positive  due  to  a  lower  frequency  component  (that  of 
seasonal  variation,  see  Figure  3.10).  Thus,  a  time  shift  of  3 
months  or  greater  would  be  needed  to  show  the  seasonal  variation 
within  the  autocorrelation.  However,  disregarding  the  positive 
nature  of  these  values,  it  is  seen  that  the  main  component  ap- 
pearing is  that  of  the  24  hour  diurnal  nature  of  Kansas  load 
demand.  In  addition,  the  weekly  demand  cycle  can  also  be  seen  as 
the  curve  maximums  dip  during  midweek  time  shift,  and  rise  for 
zero  and  one  week  time   shifts. 

Next,  the  wind  power  autocorrelation  for  both  Tuttle  Creek 
and  Plainville  shows  more  of  the  random  result  that  would  be 
expected  of  power  produced  from  wind.  As  seen  in  the  plots  for 
these  two  sites,  both  have  the  dramatic  peak  which  occurs  at  zero 
time    shift   for    autocorrelations   of   a  random   nature.      However,     it 
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is  also  seen  that  wind  is  not  completely  random  (gradual  drop 
from  100*  correlation),  and  has  80%  correlation  at  time  shifts  of 
+  2  hours.  Thus,  the  autocorrelation  shows  that  variation  in 
winds  must  occur  over  higher  time  periods  since  correlation  is 
still  around  20%  for  even  a  daily  shift.  This  is  valuable  infor- 
mation, and  is  discussed  in  Data  Validation.  Finally,  the  four 
to  five  day  frontal  cycle  is  evident  in  the  lower  frequency 
component  that  begins  to  show  two  days  from  zero  time  shift. 
Notice  that  between  Tuttle  Creek  and  Plainville,  the  frequency  of 
this  component  is  relatively  the  same,  but  shows  to  have  very 
dissimilar  waveforms. 
Crosscorrelation 

Crosscorrelations  are  computed  for  all  sites,  and  help  to 
show  how  a  particular  site  correlates  with  the  Kansas  load  de- 
mand. Note  that  all  of  the  results  from  this  calculation  have 
waveforms  which  consist  of  very  small  correlation  values.  This 
should  be  the  case  since  wind  as  noted  before  is  a  relatively 
random  signal  and  will  have  low  correlation  with  the  demand  of 
Kansas  people  and  industry.  However,  of  more  importance  is  the 
shape  that  the  correlation  waveform  takes  as  time  shift  is  in- 
creased. 

Beginning  with  Tuttle  Creek,  the  crosscorrelation  shows  that 
both  Tuttle  Creek  and  the  utility  have  a  good  24  hour  correlation 
component.  Also,  if  the  wind  profile  were  shifted  back  in  time 
by  three  hours  (later  wind  speeds  occurring  three  hours  earlier), 
then   the   24  hour    componenet   would  overlap   exactly,    and  cause    the 


peak   close    to    the    center    of    the   plot    to    indeed   occur    at    the 
center. 

Plainville  on  the  other  hand,  has  a  greater  time  shift  till 
overlap  of  the  diurnal  cycles.  The  plot  for  this  crosscorrela- 
tion  shows  Plainville  needing  a  six  hour  time  shift  to  align  the 
24  hour  cycle  of  the  wind  with  that  of  the  utility.  Note  that  in 
this  plot,  the  24  hour  cycle  is  not  as  smooth  as  that  for  Tuttle 
Creek,  and  would  thus  indicate  that  Plainville  does  not  follow 
the   diurnal    cycle   as    readily  as  Tuttle  Creek. 

Moving  to  Atlanta,  its  crosscorrelation  with  the  utility 
appears  to  be  similar  to  that  of  Plainville  (not  as  smooth  a 
diurnal  correlation  as  Tuttle  Creek),  but  shows  a  much  greater 
time  shift  of  12  hours  until  the  24  hour  cycle  has  a  peak  in 
correlation.  This  tends  to  indicate  that  Atlanta  wind  speeds 
vary  180  degrees  out  of  phase  with  the  utility,  and  have  some 
form   of   a    nocturnal    jet. 

Completing  the  single  site  crosscorrelations,  Wright  like 
Tuttle  Creek  shows  close  correlation  with  the  utility  24  hour 
cycle.  However,  it  is  seen  that  at  the  peak  of  this  correlation, 
Wright  takes  a  slight  dip  in  the  correlation  waveform  for  all  24 
multiples  of  time  shift.  This  could  be  a  characteristic  of 
Wright  itself,  or  as  with  Atlanta,  could  be  a  circumstance  of  the 
incomplete  data  set  used  in  this  analysis.  Whatever  the  case, 
Wright's  crosscorrelation  plot  shows  decreasing  wind  power  when 
the   utility    is   at    its   peak  power  production  for   the   day. 

Finally,  to  complete  the  crosscorrelation  analysis,  multiple 
site   wind  arrays   are   used   to   correlate   with    the   utility.     As    seen 
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for  a  two  site  wind  array  (Tuttle  Creek  and  Plainville  50  m),  the 
24  hour  correlation  cycles  start  to  fade,  but  show  a  four  hour 
time  shift  from  zero  to  peak  the  crosscorrelation.  For  the  three 
site  wind  array  (Tuttle  Creek,  Wright,  and  Plainville  30  m),  the 
24  hour  cycle  fades  even  more  and  now  is  hardly  recognizable. 
From  this,  it  is  noted  that  increasing  the  number  of  sites  de- 
creases the  diurnal  correlation  between  utility  and  wind  array, 
thus  causing  the  data  sets  to  tend  toward  independence  in  this 
area. 


CHAPTER  V 

CONCLUSIONS  AND  RECOMMENDATIONS 

From  the  results  and  discussion  of  ramping  effects  and 
correlation  analysis  given  in  the  previous  chapter,  the  following 
conclusions  and  recommendations  are  made. 

1.  All  sites  used  in  this  study  show  detrimental  effects  to  the 
utility  in  some  way,  shape,  or  form.  However,  at  the  50  m  level, 
Tuttle  Creek  and  Plainville  prove  to  be  better  sites  since  Atlan- 
ta showed  to  have  the  most  severe  worst  case  ramping  conditions 
of  the  three  sites.  At  30  m,  all  sites  behave  relatively  the 
same,  thus  giving  no  best  case  condition. 

2.  Since  results  for  all  sites  showed  lowering  of  the  ramping 
distributions  for  increases  in  u  and  u  ,  for  the  options  used, 
uc  =  7.0  mps  and  ur  =  14.0  mps  would  be  recommended.  Even  though 
this  selection  decreases  total  wind  energy  production,  its  main 
advantage  is  in  decreasing  excessive  ramping  due  to  low  wind 
speed  cut  in. 

3.  In  recaping  previous  results,  it  was  noted  that  as  different 
sites  were  added  to  the  utility,  that  each  site  created  a 
modified  version  of  the  utility's  average  24  honr  day.  However, 
as  seen  in  later  results,  as  uc  and  u  ,  penetration,  and  ramping 
control  were  varied,  no  noticeable  trends  varied  from  these 
modified  versions.  This  would  support  the  conclusion  that  the 
mean  24  hour  averages  are  functions  of  site  and  not  of  parameters 
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used   to  vary   the   means   of   wind  power  production. 

4.  As  expected,  increasing  penetration  causes  a  proportional 
increase  in  ramping  deviations  and  worst  case  minimum*  and  maxi- 
mums. However,  as  seen  in  the  aggregate  control  example,  a  wind 
array  of  three  sites  with  15%  penetration  can  have  relatively  the 
same  impact  as  that  of  a  5%  single  site  case,  if  control  parame- 
ters   are    chosen   correctly. 

5.  Since  increasing  ramping  control  causes  the  affected  negative 
utility  distribution  to  come  closer  to  the  utility  norm,  ramping 
control  is  recommended.  Even  though  this  might  consist  of  a 
complex  network,  such  control  would  benefit  utility  operation.  A 
typical  value  for  this  thesis  which  gives  good  control  and  mini- 
mal power  loss    is   1.0    MW/min. 

6.  In  terms  of  the  number  of  sites  in  a  wind  array,  two  or  more 
sites  per  array  is  suggested  to  smooth  the  effects  of  injected 
wind  power.  Even  though  the  aggregate  wind  array  showed  results 
of  producing  an  increase  in  the  evening  pickup  for  24  hour  means, 
this  is  of  little  importance  as  compared  to  the  resultant  de- 
crease in  worst  case  minimums  and  maximums  to  the  level  of  the 
utility  norm.  Thus,  multiple  site  farms  do  change  the  average 
operation  of  the  utility,  but  offer  large  benefits  in  terms  of 
worst   case    ramping. 

7.  For  the  correlation,  it  is  noted  that  wind  power  and  Kansas 
load  demand  do  show  a  correlative  component  on  the  diurnal  cycle. 
Even   though   the  correlation  values   are   small    this   correlation 
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pattern  is  still  very  evident.  Tuttle  Creek  and  Plainville  have 
the  closest  correlation  with  Kansas  load  demand,  while  Atlanta 
tends    to   be   12   hours   out    of    correlation. 

8.  Finally,  as  more  sites  are  added  into  a  wind  array,  the  24 
hour  correlation  cycle  decreases,  and  as  seen  in  a  three  site 
correlation  with  the  utility,  is  hardly  recognizable.  This  re- 
sult shows  that  wind  power  and  Kansas  load  demand  become  increas- 
ingly independent  in  terms  of  diurnal  correlation  as  the  number 
of    sites    increases. 

From  these  conclusions  and  recommendations,  it  can  be  seen 
that  the  negative  effect  of  added  wind  generation  can  be  reduced 
to  controllable  levels.  To  completely  see  the  benefit  of  wind 
power  however,  these  effects  would  have  to  be  incorporated  into 
a  study  on  the  benefits  and  cost  of  such  wind  power  production, 
which  lies  beyond  the  scope  of  this  thesis.  Finally,  to  verify 
the  results  for  Atlanta  50  m,  a  follow  up  study  is  recommended 
for  the  years  1981  to  1983,  when  adequate  utility  data  has  become 
available. 
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c 

C  COMBINED  RAMPING  RATE  CALCULATION  ROUTINE 

C 

C  DG  FORTRAN  5  SOURCE  FILENAME:                         RAMP3.FR 

C 

C  DEPARTMENT  OF  ELECTRICAL  ENGINEERING    KANSAS  STATE  UNIVERSITY 

C 

C  REVISION       DATE                 PROGRAMMER 

C      

c  00.0                         JULY  10,  1984                       CHRISTOPHER  DUFFEY 
C 

C 

C  PURPOSE 

C 

C  THIS  ROUTINE  CALCULATES  RAMPING  RATES  FOR  UTILITY  DATA 

C  THAT  HAS  BEEN  INJECTED  WITH  WIND  POWER, 

C  CALCULATED  THROUGH  A  WIND  POWER  PRODUCTION  MODEL. 

C 

C  ROUTINE (S)    CALLED  BY  THIS  ROUTINE 

C 

C  VMEM  -  DG  SYSTEM  CALL 

C  MAPDF  -  DG  SYSTEM  CALL 

C  CHECK  -  DG  SYSTEM  CALL 

C  VDUMP  -  DG  SYSTEM  CALL 

C  VLCAD  -  DG  SYSTEM  CALL 

C  ISTASH  -  STORE  EXTENDED  INTEGER  VALUE  ROUTINE 

C  IEXTD  -  RECALL  EXTENDED  INTEGER  VALUE  ROUTINE 

C  RINPUT  -  REAL  VARIABLE  PROMPT  ROUTINE 

C  IINPUT  -  INTEGER  VARIABLE  PROMPT  ROUTINE 

C  HNAME  -  HFTLE  FILENAME  GENERATION  ROUTINE 

C  STAT  -  DG  SYSTEM  CALL 

C  OPEN  -  DG  SYSTEM  CALL 

C  READRW  -  DG  SYSTEM  CALL 

C  CLOSE  -  DG  SYSTEM  CALL 

C  REPLY  -  UTILITY  LIBRARY  CALL 

C  OPENW  -  UTILITY  LIBRARY  CALL 

C  OPENR  -  UTILITY  LIBRARY  CALL 

C  RESET  -  DG  SYSTEM  CALL 

C 

c 

C  NOTE  1:  This  routine  uses  extended  memory  to  perform  calculations. 

C  For  more  information,   refer  to  the  subroutines  ISTASH  and 

C  IEXID.     Also,    refer  to  the  DATA  GENERAL  FORTRAN  5  PROGRAMMER'S 

C  GUIDE,   RDOS,   starting  on  page  II-21-1. 

C 

C  NOTE  2:  To  run  this  routine  properly,  the  NOVA  computer  must  be 

C  partitioned  in  the  following  manner. 

C 

C  Execute  these  commands  from  the  background, 

c  while  in  the  RDOS  directory  DPOF. 

c  CNTRL-F 

c  SMEM  7   (CH) 

C  EXPG/E  CLI    (CR) 


c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
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NOTE  3:  To  patch  this  routine  for  operation  on  different  consecutive 
years  of  demand  data,  use  the  patch  variables  IYBEG  and  IYEND 
These  variables  represent  the  years  of  data  to  be  analyzed, 
and  are  represented  as  follows. 

1  =  1980 

2  =  1981 

3  =  1982 

4  =  1983 


NOTE  4:  Noting  that  for  days  of  the  week  Sat=l,  Sun=2,  .. 
the  start  days  of  the  years  1980,  81,  82,  83  are 
respectively  4,  6,  7,  1. 


Fri=7, 


C***  COMMON  DECLARATION  USED  IN  EXTENDED  MEMORY  ALLOCATION 

COMMON.   WINDOW  (1024) 
C***  VARIABLE  DECLARATIONS 


LOGICAL       ANS 

LOGICAL       DEBUG 

LOGICAL       ANStfER 

LOGICAL       NODAT 

LOGICAL       WSKIP 

LOGICAL       ROTTRL 

REAL       SUMR(24)/24*0.0/ 

REAL       SUMRSQ(24)/24*0.0/ 

REAL       RMAX(24)/24*-999.0/ 

REAL       RKEN(24)/24*999.0/ 

REAL       STNDV(24)/24*0.0/ 

REAL       MSAN(24)/24*0.0/ 

REAL      DATA(744)/744*0.0/ 

REAL       INFO (744) /744*0.0/ 

REAL      RCNT(24)/24*0.0/ 

REAL       RLOOP 

REAL       RLOC 

REAL       RAMP 

REAL       RDAY 

REAL       WEND  ,   WDAY 

REAL       RWEEK 

REAL       RHCCR 

REAL       SQSUMR 

REAL       SIZE 

REAL       PWR 

REAL       DC 

REAL       UR 

REAL       RTCWER 

REAL       RLOCAT 

REAL       FRAC 

REAL       RLCOP1 

REAL       CPWR 

REAL       RATED 

REAL       WPLOSS 


;USER  REPLY 

;RCOTINE  DEBUGGER  FLAG 

;USER  REPLY 

;DATA  EXISTANCE  FLAG 

;SKIP  WIND  DATA  FLAG 

;RAMP  CONTROL  FLAG 

;SUM  OF  RAMPING  RATES 

;SUM  OF  SQUARED  RAMPING  RATES 

;MAX  RAMPING  RATES 

;MIN  RAMPING  RATES 

;STND  DEV  OF  RAMPING  RATES 

;MEAN  OF  RAMPING  RATES 

;DATA  FOR  1  MONTH,  LOC,  TOWER 

; INFORMATION  HEADER 

;REAL  COUNTER  FOR  CALCULATIONS 

;REAL  LOOP  PARAMETER 

;SEAL  MEM  LOCATION  FOR  TEST 

;RAMP  RATE  FOR  1  HR  (MW/MIN) 

;NUMBER  OF  CALCULATED  DAYS 

;WEEKEND,  WEEKDAY  HOUR  TOTALS 

;NUMBER  OF  WEEKS  IN  PERIOD 

;NUMBER  OF  CALCULATED  HOURS 

,-SQUARE  OF  RAMPING  SUMMATION 

;RETURNED  SIZE  IN  BYTES 

;  CALCULATED  WINDPCWER 

;  CUT- IN  WINDSPEED 

,-RATED  WIND  SPEED 

;50  OR  30  METER  DATA 

;SITE  LOCATION  NUMBER 

;FRACTION  OF  WIND  POWER 

;REAL  LOOP  PARAMETER 

;  COMBINED  POWER  EFFECTS 

;RATED  WINDFARM  OUTPUT  (MW) 

;SUM  OF  LOST  WIND  POWER 


REM, 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER      WEEK 

INTEGER       TWK 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 


CUT 

WPSUM 

MAXRAMP 

OJTINWS(8)/8*0.0/ 

RATEDWS(8)/8*0.0/ 

NSUM(51)/51*0.0/ 

PSUM(51)/51*0.0/ 

NRES(  51)/51*0.0/ 

PRES(  511/51*0.0/ 

LLIMTT 

ULIMTT 

PCNT 

NOT 

TRAMP 

RATE1    ,   RATE2 

FRACTIDN(8)/8*0.0/ 

MNUM(4, 121/48*0/ 

IYBEG  ,    IYEND 

NUMDAY 

IDAY  ,    IMONTH 

IHCXJR  ,    IYEAR 

PWR1    ,    PSR2 

PDIFF 

N(12)/12*0/ 

SDOY(4)/4, 6,7,1/ 


BCNT 

IEEC 

IPWR 

ICHK(18) 

IPARM 

NPWR2 

NUMSITE 

SITE  (81/8*0/ 

TCWER( 8)/8*0/ 

I  ,   J 

LQAD1   ,   LQAD2 
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;  COUNTER  FOR  GOOD  DATA 

;TCTAL  WIND  POWER  PRCD 

,-MAXIMUM  CONTROL  RAMPING 

;CUT-IN  WINDSPEEDS 

; RATED  WINDSPEEDS 

;NEG  RAMPING  HIST  SUMS 

;POS  RAMPING  HIST  SUMS 

;NEG  DISTRIBUTION  %  RESULTS 

;POS  DISTRIBUTION  %  RESULTS 

;LCWER  BIN  LIMIT  (HIST) 

;UPPER  BIN  LIMIT  (HIST) 

;#  OF  PCS  RAMP  RATES 

;#  OF  NEG  RAMP  RATES 

;TEMP  -RAMP 

;PRINT  OUT  RAMP  RATES 

.-FRACTION  OF  FARM  POWER 

;DAYS  QJ  EACH  MONTH 

;YEAR  LOOP  PARAMETERS 

;NUMBER  OF  TOTAL  CRUNCH  DAYS 

;LOOP  PARAMETERS 

;LCOP  PARAMETERS 

;TWO  CCNSEC  HOUR  POWERS  X100 

;ENERGY  DIFFERENCE  (MW  X100) 

;TOTAL  HOURS/MONTH  FOR  3  YEARS 

;START  DAY  OF  YEAR  80,81,82,83 

;NUMBER  OF  WEEKS  IN  PERIOD 

;LCOP  VARIABLE  FOR  WEEK 

;RETURNED  BLOCK  COUNT 

;RECORD  NUMBER  FOR  DATA 

;INTEG  VAL  OF  WIND  POWER  X100 

;STATUS  VARIABLE 

;LCOP  PARAMETER 

;NEW  PWR2  FROM  RAMP  CONTROL 

,-NUMBER  OF  SITES 

;SITE  NUMBERS 

; TOWER  HEIGHTS 

;LCOP  PARAMETERS 

;TWO  CONSECUnVE  LOADS 

;H-FILE  DATA  FILENAME 


DOUBLE  PRECISION  COMPLEX       FILNAM 

C***  VARIABLE  INITIALIZATIONS 

IYBEG     =  1 
IYEND     »  3 
DEBUG     =   .FALSE. 
ANSWER  =   .FALSE. 

C***  nnnALIZE  EXTENDED  MEMORY  FOR  REAL   (4-BYTE)    VARIABLE  SIZE 

CALL  VMEM      (ICNT  ,    IERR) 

CALL  MAPDF   (ICNT+2   ,   WINDOW   ,   2   ,    IERR) 

CALL  CHECK  (IERR) 

CALL  MAPDF  (2  ,  IERR) 

CALL  CHECK  (IERR) 

C***  INITIALIZE  DAYS  IN  MONTH  ARRAY   (FOR  I,   1=80,2=81,3-82,4=83) 


DO  10  I  =  1  ,  4  91 

MNUM(I,1)  =  31 
HNUM(I,2)  =  28 
MNUM(I,3)  =  31 
MNUM(I,4)  =  30 
MNUM(I,5)  =  31 
MNDM(I,6)  =  30 
MNUM(I,7)  =  31 
MNDM(I,8)  =  31 
MNUM(I,9)  =  30 
MNUM(I,10)  =  31 
MNUM(I,11)  =  30 
MNUM(I,12)  =  31 
10     CONTINUE 

C***  ACCOUNT  FOR  LEAP  YEAR  (1980) 

MNUM(1,2)  =  29 

C***  ASK  USER  TOR  DEBUGGING  OPTION  AND  EXT  MEM  LOAD 

CALL  REPLY   ("DEBUG  ?    (YES/NO)    ■   ,    DEBUG) 
CALL  REPLY   ("SKIP  AVE  CALCULATIONS  ?    (YES/NO)    "   ,   WSKLP) 
CALL  REPLY    ("LOAD  EXT  MEM  FROM  DISK  ?    (YES/NO)    ■       ANSWER) 
IF   (.IDT.   ANSWER)    GO  TO  14  .  Uts/Nu)        '   ™SWER) 

C***  LOAD  IN  AN  EXTENDED  MEMORY  FILE 

CALL  OPENR   (3  ,  "INPUT  FILE  ?   "   ,   0    ,   SIZE) 

CALL  VLOAD    (3  ,  BCNT  ,    IERR) 

CALL  CLOSE   (3  ,  IERR) 

14  CONTINUE 

C***  INITIALIZE  THE  NUMBER  OF  DAYS  TO  BE  USED  IN  CALCULATIONS 

NUMDAY  =  0 

CO  17  I  =  IYBEG  ,  IYEND 
DO  15  J  =  1  ,  12 

NUMDAY  =  NUMDAY  +  MNUM(I,J) 

15  CONTINUE 
17      CONTINUE 

C***  INITIALIZE  THE  2ND  HALF  OF  EXTENDED  MEMORY  WITH  0 'S 
TYPE   "   " 

TYPE   "INITIALIZING  2ND  HALF  OF  EXTENDED  MEMORY" 

RLOOP  =  26500.0 

CO  13  I  »  1   ,   NUMDAY  *  24  +  1 

RLOOP  =  RLOOP  +  1.0 

CALL  ISTASH   (0   ,   RLOOP) 
13  CONTINUE 

Zll  ^™  P  WIND  mTA  POR  THE  GIVEN  SITE,   AND 
C***  STORE  THE  CALCULATED  POWER  IN  EXTENDED  MEMORY 

C***  PROMPT  THE  USER  FOR  INFORMATION 


NUMSITE  =0  92 

THE   "    " 

CALL  ROBOT   (RATED  ,   15    ,    "RATED  WINDFARM  POWER  COTPOT  1  "1 

^.^LS^tcTdg  0DmK1  ?  (YES/N0> "  - m) 

CALL  RINFUT   (MAXRAMP  ,    14    ,    "MAX  RAMPING  RATE  ?    (MS/MIN)    ») 

16  CALL  RINHJT  (RLOCAT  ,    6   ,    "LOCATION  ?     ") 

CALL  RINFOT  (RTCWER  ,   9   ,    "30  OR  50  METER  ?  ") 

21  CALL  RINFDT  (DC  ,   10   ,    "CDT-IN  WIND  SPEED  ?  ") 

CALL  RDJPDT  (DR  ,   10    ,    "RATED  WIND  SPEED  I  ") 

CALL  RINFDT  (FRAC  ,   19    ,    "FRACTION  OF  TOTAL  POWER  PRODUCTION  ?     ") 

NOMSITE  =  NDMSITE  +  1 
SITE (NUMSITE)    =  IFIX(RLOCAT) 
TOWER (NDMSITE)    =  IFIX (RTCWER) 
CUTCNWS  (NDMSITE)    =  UC 
RATEDWS(NUMSITE)    =  UR 
FRACTION(NDMSITE)    =  FRAC 

RLCOP  =  26500.0 
NODAT  =   .TRUE. 

DO  50   HEAR  =  IXBEG  ,    IYEND 
DO  40   IM3NTH  =  1   ,  12 

CALL  HNAME    (IMONTH   ,    HEAR  ,   FILNAM) 

C***  CHECK  FOR  DEBUG  -  FILENAME  CHECK 

IF(.NCT.   DEBUG)    GO  TO  19 
WRITE (10, 18)    FILNAM 

18  FORMAT(A16) 

19  CONTINUE 

C***  CHECK  FOR  FILE  EXISTANCE 

CALL  STAT   (FILNAM  ,    ICHK  ,    IERR) 
IF   (IERR  .£Q.   13)    GO  TO  30 

C***  THIS  FILE  EXISTS,   READ  IN  1HE  INFORMATION  HEADER  FOR  THE  FILE 

CALL  OFEN   (3   ,   FILNAM  ,   2    ,   2976   ,    IERR) 
CALL  READRW   (3,1,    INFO  ,   1   ,    IERR) 

C***  SEARCH  FOR  DATA 

I  -  1 
IREC  =  2 

20  CONTINUE 

IF(INFO(I)    .GT.  90.0)   GO  TO  28 

IFUINFO(I).NE.RLCCAT)    .OR.    (INFO(I+l)  .NE.RTCWER) )    GO  TO  25 

C***  DATA  HAS  BEEN  PCOND  FOR  THIS  HMTH,   22AR,    LOCATION,  AND  TOWER 
THE   "DATA  FODND  FOR"   ,    WMB   ,    HEAR+79 


IF    (DEBUG)    TOPE   "START  LOCATION  =   "    ,    RLOOP+1.0 
NODAT  =    .FALSE. 

CALL  READRW    (3    ,    IREC  ,    DATA  ,    1   ,    D5RR) 
GO  TO  35 
25  CONTINUE 

IREC  =  IREC  +  1 
I  »  I  +  2 
GO  TO  20 

C***  NO  DATA  FOUND 

28  CONTINUE 

TYPE   "NO  DATA  FOR  "   ,    IMDNTH   ,    IYEAR+79 
NODAT  =    .TRUE. 
GO  TO  35 

C***  NO  FILE  FOUND 

30  CONTINUE 

THE   "NO  FILE  FOR  "   ,    MONTH   ,    IYEAR+79 
NODAT  =    .TRUE. 

C***  CALCULATE  THE  ELECTRIC  WIND  POWER  AND  STORE  IN  EXTENDED  MEMORY 

35  CONTINUE 

DO  38  IHOUR  =  1   ,    MNUM(IYEAR,IMDNTH)*24 
IF(DATA(IHCUR)    .LT.   UC)    PWR  =  0.0 
IF(DATA(IHCUR)    .GT.   UR)    PWR  =  RATED  *  FRAC 
IF((DATA(IHCUR).GT.UR)    .OR.    ( DATA ( IHOUR)  .LT.UC})    GO  TO  37 
PWR  =  RATED  /   (UR  -  UC)    *   (DATA(IHCUR)    -  UC)    *  FRAC 

37  CONTINUE 

IF(DATA(IHCUR)    .GT.  90.0)   PWR  =  3200.0 
IF (NODAT)    PWR  =  3200.0 
RLOOP  =  RLOOP  +  1.0 

C***  THE  POWER  TIMES  10   IS  STORED  TO  MAINTAIN  ACCURACY 

IPWR  =  IFTX(PWR*10.0) 

IF   (IPWR  .GT.   30000)   GO  TO  39 

IPWR  =  IPWR  +  IEXTD(RLOOP) 

39  CALL   ISTASH   (IPWR  ,    RLOOP) 

38  CONTINUE 

CALL  CLOSE    (3   ,    IERR) 

40  CONTINUE 
50              CONTINUE 

CALL  REPLY   ("ADD  ANOTHER  LOCATION  ?   (YES/NO)    "   ,   ANSI 
IF   (ANS)    GO  TO  16 

C***  CHECK  FOR  DEBUG  -  MEMORY  TEST 

55  TF   (.NOT.DESUG)    GO  TO  78 

CALL  REPLY   ("MEM  TEST  ?   (YES/NO)    "    ,   ANSWER) 

IF(.NOT.ANSWER)    GO  TO  78 
60  CALL  RLNHJT   (RLOC  ,    12   ,    "REAL  LOC  FOR  MEM  TEST  ?  ") 

IGET  =  IEXID   (RLOC) 

WRITE   (10,70)    RLOC  ,    IGET 


70  FORMAT   ("MEM  LOC  (",F7.0,")    HAS  THE  VALUE  >   ",I7)  94 

C***  PROMPT  USER  FOR  CHANGE  OF  VALUE 

CALL  REPLY   ("CHANGE  THIS  MEM  LOC  ?    (YES/NO)    "    ,   ANSWER) 
IF   (.NOT.   ANSWER)    GO  TO  75 

CALL  IINFUT   (IGET  ,   10   ,    "THE  NEW  VALUE  IS  ?      ") 
CALL  ISTASH    (IGET  ,    RLOC) 
75  CONTINUE 

CALL  REPLY    ("RON  AGAIN  ?   (YES/NO)    "    ,   ANSWER) 
IF   (ANSWER)   GO  TO  60 

C***  FRCMFT  USER  TO  EUMP  EXTENDED  MEMORY  TO  DISK 

78  CALL  REPLY   ("DUMP  EXT  MEM  TO  DISK  ?    (YES/NO)    "    ,   ANSWER) 

IF   (.NOT.  ANSWER)    GO  TO  80 

CALL  OPENW    (3   ,    "OUTPUT  FILENAME  FOR  MEM  DUMP  ?   "   ,    0    ,   SIZE) 

CALL  VDUMP   (3   ,    BCNT  ,    IERR) 

CALL  CLOSE    (3   ,    IERR) 
80  CONTINUE 


C***  UPWARDS  RAMPING  CONTROL  ALGORITHM 
IP   (.NOT.   RCNTRL)    GO  TO  53 


TYPE   "   " 

TYPE   "PERFORMING  RAMP  CONTROL" 

WPLOSS  =  0.0 

CUT  =  0.0 

RLOOP  =  26500.0 

WFSUM  =0.0 

IF   dEXTD(26501.0) .LT.30000)   WPSUM  =  FLOAT (IEXTD   (26501.0)1/10.0 

DO  52  IPARM  =  1   ,   NUMDAY*24 
RLOOP  =  RLOOP  +1.0 
PWR1  =  IEXTD   (RLOOP) 
PWR2  =  IEXTD   (RLOOP  +  1.0) 

IF    (PWR2 .LT.30000)   WPSUM  =  WPSUM  +  FLOAT(PWR2)/10  0 
IF   (PWR2.LT.30000)   CNT  =  GST  +  1.0 
IF  ((PWR1.GT.30000).OR. (PWR2.GT.30000) )   GO  TO  52 
PDIFF  =  EWR2  -  IWR1 
RAMP  =  FLOAT(PDIFF)/600.0 
IF   (RAMP.LE.MAXRAMP)    GO  TO  52 
NIWR2  =  PWR1  +  IFIX(MAXRAMP*600.0) 
CALL  ISTASH    (NPWR2   ,   RLCOP+1.0) 
WPLOSS  =  WPLOSS  +  FLOAT(PWR2-NPWR2)/10.0 
IF   ((PWR2-NIWR2).LT.0)    TYPE   "RAMP  CONTROL  LOGIC  ERROR" 

52  CONTINUE 

53  OONTTNUE 

C***  CALCULATE  COMBINED  POWER  PRODUCTION 

TYPE  "  " 

TYPE  "COMBINING  WHS  UTILITY  POWER  DATA" 

RLOOP1  =  0.0 

RLOOP    =  26500.0 


95 

DO  84   IYEAR  =   IYBEG   ,    IYEND 
DO  83  IMONTH  =  1   ,   12 

DO  82   HAY  =  1   ,    MNUM(IYEAR  ,    IMONTH) 
DO  81   IHOUR  =  1   ,   24 

RLOOP1  =  ELOOP1  +1.0 
SLOOP     =  RLOOP     +  1.0 

CPWR  =  nXWr(IEXm(PLOOPl))-FLQAT(IEXID(ELCOP))/10.0 
IF   (IEXTD(RLOOP)    .CT.   30000)    CPWR  =  32000.0 
CALL  ISTASH   (IFIX(OWR)    ,   RLOOP1) 

81  CONTINUE 

82  CONTINUE 

83  CONTINUE 

84  CONTINUE 

C***  SET  LAST  ELEMENT+1  ID  32000 

CALL  ISTASH   (32000   ,   RLOOP1  +  1.0) 

C***  PROMPT  USER  TO  DUMP  EXTENDED  MEMORY  TO  DISK- 
CALL  REPLY   ("DUMP  EXT  MEM  TO  DISK   (YES/NO)    ?   "   ,   ANSWER) 
IF   (.NOT.   ANSWER)    GO  TO  86 

CALL  OPEtW   (3   ,    "OUTPUT  FILENAME  FOR  MEM  DUMP  >   "   ,   0    ,    SIZE) 

CALL  VDUMP   (3   ,    BCNT  ,    IERR) 
CALL  CLOSE   (3   ,    IERR) 

86  CONTINUE 

C**»  FIND  THE  DISTRIBUTION  OF  THE  RAMPING  RATES 

TYPE   "    " 

TYPE   "PERFORMING  DISTRIBUTION  CALCULATIONS" 

RLOOP1  =  26500.0 

RLOOP  =0.0 

POST  =  NQJT  =  0.0 

DO  92  I  =  1  ,  50 

NSUM(I)    =0.0 

PSUM(I)    =0.0 

NRES(I)    =  0.0 

PRES(I)   =  0.0 
92  CONTINUE 

DO  95   IPARM  =  1   ,   NUMDAY*24 
RLOOP  =  RLOOP  +  1.0 
LOAD1  =  IEXH)   (RLOOP) 
LOAD2  =  lEXTD   (RLOOP  +  1.0) 

IF  ((LOAD1.GT.30000).OR.(LQAD2.GT.30000))   GO  TO  95 
HJIFF  =  LQAD2  -  LOAD1 
RAMP  =  FLOAT(PDIFF)/60.0 

C*»*  GENERATE  THE  HISTOGRAM  OF  THE  DATA  SET 

DO  94  I  =  1   ,   51 

IF  (RAMP  .LE.  -0.10)   GO  TO  93 
LLDOT  =  FLOAT(I-1)/5.0  -  0.10 
ULIKtT  =  FLOAT(I-1)/5.0  +  0.10 
IF((RAMP.LE.LLIMIT).CR.  (RAMP.GT.ULIMIT) )    GO  TO  94 


96 

PCNT  =  POTT  +  1.0 
PSUM(I)  =  PSUM(I)  +  1.0 
IF  (RAMP  .LE.  0.10)  GO  ID  93 
GO  TO  95 

93  CONTINUE 
TRAMP  =  -RAMP 

LLDCCT  =  FLOAT(I-1)/5.0  -  0.10 
ULIMIT  =  FLOAT(I-1)/5.0  +  0.10 

IF((TRAMP.LT.LLIMTT).aR. (TRAMP. GE. ULIMIT) )    GO  TO  94 
NCNT  =  NOT  +  1.0 
NSUM(I)    =  NSUM(I)    +  1.0 
GO  TO  95 

94  CONTINUE 

95  CONTINUE 

C***  CHECK  FOR  DEBUS  -  NSUM  AND  PSUM 

IF  (DEB03)  TYPE  "NSDM  AND  PSUM" 
IF  (DEBUS)  TYPE  NSUM 
IF  (DEB03)  TYPE  PSUM 

C***  SUM  THE  HISTOGRAM  INTO  A  DISTRIBUTION 

DO  97  I  -  1  ,  51 

DO  96  J  =  52-1  ,  51 

NRES(52-I)  =  NRES(52-I)  +  NSUM(J) 
PRES(52-I)  =  PRES(52-I)  +  PSUM(J) 

96  CONTINUE 

97  CONTINUE 

C***  TURN  THE  SUMS  INTO  PERCENT  OF  THE  TIME  THAT  THIS 
C***  RAMPING  RATE  OR  GREATER  HAS  BEEN  OBSERVED 

DO  98  I  -  1  ,  51 

NRES(I)  =  (NRES(I)/NCNT)*100.0 
PRES(I)  =  (PRES(I)/PCNT)*100.0 

98  CONTINUE 

C***  PERFORM  CALCULATIONS  FOR  24  HOUR  AVERAGES  OVER  DESIRED  PERIOD, 
C***  CALCULATE  MEANS,  STANDARD  DEVIATIONS,  MTNS,  MAXS 

IF  (WSKIP)  GO  TO  450 
RLOOP  =0.0 
DO  85  I  -  1  ,  24 
RONT(I)  =  0.0 
85     CONTINUE 

TYPE  "  " 

TYPE   "PERFORMING  24  HOUR  COMPUTATIONS" 
DO  100   HAY  =  1   ,   NUMDAY 
DO  90  IHCUR  »  1   ,  24 

RLOOP  =  RLOOP  +  1.0 

PSF1  »  rEXTD   (RLOOP) 

PHR2  »  IEXTD   (RLOOP  +  1.0) 

IF((PWR1.CT.30000)    .OR.    (PWR2.GT.30000) )   GO  TO  90 

PDIFF  =  EWR2  -  PSR1 


RAMP  =  FLOAT(PDIFF)    /  60.0 

RMAX  (IHCOR)    =  AMAX1    (RMAX(IHOUR)    ,   RAMP) 

RMINdHCOR)    =  Alma    (RMTN(IHOUR)    ,   RAMP) 

SUMR(IHOUR)    =  SUMR( IHCOR)    +  RAMP 

SUMRSQt  IHCOR)    =  SOMRSQ ( IHOOR)    +  RAMP*RAMP 

RCNT(IHCUR)    =  RCNT(IHOUR)    +1.0 

90  CONTINUE 

100  CONTINUE 

C***  FINISH  24  HOUR  CALCULATIONS 

DO  110  IHOUR  =  1   ,  24 
REAY  =•  RQJT(IHOUR) 
MEAN(IHCCR)    =  SUMR(  IHCOR) /REAY 
SQSOMR  =  SUMR(  IHCOR)    *  SUMR(IHCUR) 

STNDV(IHCOR)    =   (SUMRSQ(IHOUR)-SQSUMVREAY)/(RDAY-1.0) 
STNEV(IHCOR)    =  SQRT   (STNDV(  IHCOR) ) 
110  CONTINUE 

C***  CHECK  FOR  DEBUG  -  24  HCOR  TEST 

IF(.N0T.  DEBUS)  GO  TO  112 

TYPE  "MEANS  FOR  24  HOUR" 

TYPE  MEAN 

TYPE  "STNDVS  FOR  24  HOUR" 

TYPE  STNEV 

TYPE   "MDJS  FOR  24  HOUR" 

TYPE  RMIN 

TYPE  "MAXS  FOR  24  HOUR" 

TYPE  RMAX 

PAUSE  "PRESS  ANY  KEY  TO  CONTINUE" 

112  CONTINUE 

C***  WRITE  24  HCOR  DATA  TO  THE  LINE  PRINTER 

WRITE  (12,299) 

DO  113  I  =  1  ,  NUMSITE 

WRITE  (12,300)  SITE(I),TCWER(I),CUTINWS{I),RATECWS(I), 
FRACTION(I)*100.0 

113  CONTINUE 

WRITE  (12,590)  RATED 

IF  (RQJTRL)  WRITE (12,595)  MAXRAMP 

WRITE  (12,301) 

DO  114  I  =  1  ,  12 

WRITE  (12,302)  I  ,  MEAN(I)  ,  1+12  ,  MEAN(I+12) 

114  CONTINUE 

WRITE    (12,303) 
DO  116  I  =  1   ,   12 

WRITE   (12,302)    I   ,   STNDV(I)    ,    1+12   ,   STNDV(I+12) 
116  CONTINUE 

WRITE    (12,299) 

DO  113  I  -  1   ,   NUMSITE 

WRITE   (12,300)    SITE(I),rCWER(I),CUTIlWS(I),RATEDWS(I), 

no  rr™*™,,  FRACTION(I)*100.0 

119  CONTINUE 
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WRITE  (12,590)  RATED 

IF  (RCNTRL)  WRTTE(12,595)  MAXRAMP 

WRITE  (12,304) 

DO  117  I  =  1  ,  12 

WRITE  (12,302)   I  ,  RMDJ(I)  ,  1+12  ,  RMIN(I+12) 
U7    CONTINUE 

WRITE  (12,305) 
DO  118  I  -  1  , 12 

WRITE  (12,302)  I  ,  RMAX(I)  ,  1+12  ,  RMAX(I+12) 

118        continue 

299  format ("1", ////////,"  ",23x,  "combined  ramping  computations",//) 

300  forhatt"  ",9x, "location  ",12,"  ",i2,"  m,  uc=",f4.1, 
*  "  mps,  ur=",f4.1,"  mps,  %  farm  =  ",f5.1) 

301  format(////,"  ",32x,"24  hour  means",/) 

302  FORMATt"   ",18X,I3,"  — >  ",F6.3,10X,I3,"  — >   ",F6.3) 

303  FORMAT(////,"  ",24X,"24  HCOR  STANDARD  DEVIATIONS",/) 

304  FORMATt////,"  ",30X,"24  HCOR  MTNIMUMS",/) 

305  FORMATt////,"  ",30X,"24  HCOR  MAXIMUMS",/) 

C***  WRTTE  24  HCOR  DATA  TO  DISK 

TYPE  "  " 

TYPE  "WRITING  24  HOUR  DATA  TO  DISK" 

CALL  OPENW    (1   ,    "FILENAME  FOR  MEAN  24  HOURS  ?   "  ,  96    ,  SIZE) 

CALL  OPENW   (2   ,    "FILENAME  FOR  STDV  24  HOURS  ?   "  ,  96    ,  SIZE) 

CALL  OPENW    (3   ,    "FILENAME  FOR  RMIN  24  HOURS  ?   "  ,  96    ,  SIZE) 

CALL  OPENW    (4   ,    "FILENAME  FOR  RMAX  24  HOURS  ?   "  ,  96    ,  SIZE) 

CALL  WRTTRW  (1,1,  MEAN  ,   1   ,    IERR) 

CALL  WRTTRW  (2,1,  STNDV  ,   1   ,    IERR) 

CALL  WRTTRW  (3,1,  RMTN  ,   1   ,    IERR) 

CALL  WRTTRW  (4,1,  RMAX  ,    1   ,    TERR) 

CALL  CLOSE  (1  ,  IERR) 

CALL  CLOSE  (2  ,  TERR) 

CALL  CLOSE  (3  ,  TERR) 

CALL  CLOSE  (4  ,  TERR) 

C***  INITIALIZE  ARRAYS 

DO  115  I  =  1  ,   24 

SUMR(I)    =  SDMRSQ(I)    =  STNDV(I)    =  K£AN(I)    =  RQJT(I)    =0.0 
RMAX(I)    =  -999.0 
RMIN(I)   =  999.0 
115  CONTINUE 

C***  PERFORM  CALCULATIONS  FOR  12  MONTHS,   AVERAGE,   MEAN,   STANDARD 
C***  DEVIATION,    MTN  AND  MAX  FOR  EACH  MONTH. 

RLCOP  =0.0 

TYPE   "   " 

TYPE   "PERFORMING  12  MONTH  CALCULAnCNS" 
DO  150  IYEAR  =  IYBEG  ,    IYEND 
DO  140   IM3JTH  =  1  ,  12 


DO  130   HftY  =  1   ,   MNUM(IYEAR,IMONTH) 
DO  120   IHCCR  =  1   ,  24 
RLOOP  =  RLOOP  +1.0 
EWR1  =  IEXID    (RLOOP) 
PWR2  =  IEXTD    (RLOOP  +  1.0) 

IF((PWR1.GT.30000)    .OR.    (FWR2.GT.30000) )    GO  TO  120 
PDIFF  =  EWR2  -  EWR1 
RAMP  =  FLOAT (PDIFF)   /  60.0 
RMAX(IMONTH)    =  AMAX1    (RMAX(IMONTH)    ,   RAMP) 
RMDJ(IMONTH)    =  AMTNl    (RMINdMDNTfl)    ,   RAMP) 
SUMR(IMONTfl)    =  SUMR(IMONTH)    +  ABS(RAMP) 
SUMRSQ(IMONTH)    =  SUMRSQ(IMONTH)    +  RAMP*RAMP 
RCNT(IMONTH)    =  RCNT(IMONTH)    +  1.0 

120  CONTINUE 

130  CONTINUE 

140  CONTINUE 

150  CONTINUE 

C***  FINISH  12  MONTH  CALCULATIONS 

DO  190   IM0NIH  =  1   ,  12 
RBCOR  =  RCNT(IM3NTH) 
MEAN(IMONTfl)    =  SUMR(IMONTH)    /  RHCUR 
SQSUMR  =  SUMR(IMONTH)    *  SUMR(IMONTH) 

STNDV(IMONTH)    =   (SaMRSQ(IMDNTH)-SQSUMR/RHCUR)/(RHajR-1.0) 
STNDV(IMONTH)    =  SQRT(STNDV(IM0NTH)) 
190  CONTINUE 

C***  CHECK  FOR  DEBUG  -  12  MONTH  TEST 

IF(.NCT.  DEBUG)   GO  TO  195 

TYPE   "MEANS  FOR  12  MONTH" 

TYPE  MEAN 

TYPE   "STNEVS  FOR  12  MONTH" 

TYPE  STNDV 

TYPE  "MINS  FOR  12  MONTH" 

TYPE  RMIN 

TYPE   "MAXS  FOR  12  MONTH" 

TYPE  RMAX 

PAUSE  "PRESS  ANY  KEY  TO  CONTINUE" 

195  CONTINUE 

C***  WRITE  12  MONTH  DATA  TO  THE  LINE  PRINTER 

WRITE (12, 299) 

DO  180  I  =  1  ,  NUMSITE 

WRITE  (12,300)  SITE(I),TCWER(I),CUTTNWS(I),RATEEWS(I), 
FRACTION ( I) *100.0 
180     CONTINUE 

WRITE  (12,590)  RATED 

IF  (RCNTRL)  WRITE(12,595)  MAXRAMP 

WRITE  (12,309) 

DO  196  I  =  1  ,  6 

~.      WRITE  t12'302'  I  -  MEAN(I)  ,  1+6  ,  MEANU+6) 

196  CONTINUE 

WRITE   (12,306) 


197 


198 


DO  197  I  = 
WRITE 
CONTINUE 


1  ,  6 
(12,302) 


WRITE  (12,307) 
DO  198  I  =  1  ,  6 

WRITE  (12,302) 
CONTINUE 


100 


I  ,  STNEV(I)  ,  1+6  ,  SH!EV(I+6) 


RMIN(I)  ,  1+6 


RMIN(I+6) 


WRITE  (12,308) 
DO  199  I  =  1  ,  6 

WRITE  (12,302)  I 
199    CONTINUE 


RMAX(I)  ,  1+6  ,  RMAX(I+6) 


309 
306 
307 
308 


F0RMAT("  ",///,"  \25X,"12  MONTH  MEANS  (ABSOLUTE)",/) 

FORMAT(/,"  ",18X,"12  MONTH  STANDARD  DEVIATIONS  (ABSOLUTE) ",/) 

FORMAT(/,"  ",29X,"12  MONTH  MffllMJMS",/) 

FORMAT(/,"  ",29X,"12  MONTH  MAXIMUMS",/) 


C***  WRITE  12  MONTH  DATA  TO  DISK 

TYPE  "  " 

TYPE  "WRITING  12  MONTH  DATA  TO  DISK" 

CALL  OPENW  (1  ,  "FILENAME  FOR  MEAN  12  MONTH  ? 

CALL  OPENW  (2  ,  "FILENAME  FOR  STDV  12  MONTH  ? 

CALL  OPENW  (3  ,  "FILENAME  FOR  RMIN  12  MONTH  ? 

CALL  OPENW  (4  ,  "FILENAME  FOR  RMAX  12  MONTH  ? 

CALL  WRITRW  (1,1,  MEAN  ,  1  ,  IERR) 
CALL  WRITRW  (2,1,  SINDV  ,  1  ,  IERR) 
CALL  WRITRW  (3,1,  RMIN  ,  1  ,  IERR) 
CALL  WRITRW  (4,1,  RMAX  ,  1  ,  IERR) 

CALL  CLOSE  (1  ,  IERR) 

CALL  CLOSE  (2  ,  IERR) 

CALL  CLOSE  (3  ,  IERR) 

CALL  CLOSE  (4  ,  IERR) 

C***  INITIALIZE  ARRAYS 


48 
48 
48 
48 


SIZE) 
SIZE) 
SIZE) 
SIZE) 


200 


DO  200   I  =  1  ,  2 

SUMR(I)  =  SUMRSQl 
RMAX(I)  =  -999.0 
RMIN(I)   =  999.0 

CONTINUE 


I)    =  SINDV(I)    =  M5AN(I)    =  0.0 


C***  PERFORM  CALCULATIONS  FOR  WEEKDAY-WEEKEND. 

C 

C***  DISCARD  THE  FIRST  FEW  DAYS  UP  TO  THE  START  OF  A  NEW 

C***  WEEK  AND,   THE  LAST  FEW  DAYS  AT  THE  END  OF  THE  DATA 

C***  THAT  FORM  AN  INCOMPLETE  WEEK. 

RLCOP  =  0.0 

RDAY  =  FLOAT(NUMDAY  -  1) 

WEND  =  0.0 

WDAY  =  0.0 
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IF(SDOY(IYBEG)    .BQ.    1)    GO  TO  210 
RLCOP  =  8.0  -  FLOAT  (SDOY(IYBEG)) 


210  CONTINUE 


{WEEK  =   (RDAY  -  RLCOP)    /  7.0 
WEEK  =  IFIX(RWEEK) 


H      ■ 


TYPE 

TYPE   "PERFORMING  WEEKEND-WEEKDAY  CALCULATIONS" 

DO  260  IWK  =  1  ,  WEEK 

C***  PERFORM  WEEKEND  TOTALS 

DO  230  IDAY  -1,2 

DO  220   IHCUR  =  1   ,   24 

RLCOP  =  RLCOP  +  1.0 

PWR1  »  IEXTD   (RLCOP) 

EWR2  =  IEXID   (RLCOP  +  1.0) 

IF((PWR1.GT.30000)    .OR.    (PWR2.GT.30000) )    GO  TO  220 

PDIFF  ■  PWR2  -  PWR1 

RAMP  =  FLCAT(PDIFF)    /  60.0 

RMAX(l)    =  AMAXl    (RMAX(l)    ,   RAMP) 

RMIN(l)    ■  AMIN1    (RMIN(l)    ,   RAMP) 

SUMR(l)    =  SOMR(l)    +  ABS(RAMP) 

SUMRSQ(l)    =  SUMRSQ(l)    +  RAMP*RAMP 

WEND  =  WEND  +  1.0 
220  CONTINUE 

230  CONTINUE 

C***  PERFORM  WEEKDAY  TOTALS 

DO  250   IDAY  =3,7 

DO  240  IHCUR  =  1  ,  24 

RLCOP  =  RLCOP  +  1.0 

PWR1  =  IEXTD   (RLCOP) 

PSR2  =  IEXID   (RLCOP  +  1.0) 

IF((PWR1.GT.30000)    .OR.    (PWR2.GT.30000) )   GO  TO  240 

PDIFF  =  PWR2  -  PWR1 

RAMP  =  FLQAT(PDIFF)   /  60.0 

RMAX(2)    =  AMAXl    (RMAX(2)    ,   RAMP) 

RMIN(2)    =  AMDn.    (RMIN(2)    ,   RAMP) 

SUMR(2)    =  SUMR(2)    +  ABS(RAMP) 

SUMRSQU)    =  SUMRSQ(2)    +  RAMP*RAMP 

WDAY  -  WDAY  +  1.0 
240  CONTINUE 

250  CONTINUE 

260  CONTINUE 

C***  FINISH  WEEKEND-WEEKDAY  CALCULATIONS 

MEAN(l)    =  SUMR(l)    /  WEND 

MEAN(2)    =  SUMR(2)   /  WDAY 

SQSUMR  »  SUMR(l)    *  SUMR(l) 

STNDV(l)    =   (SUMRSQ(l)-SaSUMR/WEND)/(WEND-1.0) 

STNDV(l)    -  SQRT(SOCV(l)) 

SQSUMR  »  SUMR(2)    *  SUMR(2) 


STNDV(2)    =   (SUMRSQ(2)-SQSUMR/WDAY)/(WDAY-1.0) 
STNDV(2)    =  SQRT(STNDV(2)! 

C***  CHECK  FOR  DEBUG  -  WEEKEND-WEEKDAY  TEST 

IF (.NOT.   DEBUS)   GO  TO  265 

TYPE   "WEEKDAY-WEEKEND  MEANS" 
TYPE  MEAN(2)    ,   MEAN(l) 
TYPE   "WEEKDAY-WEEKEND  STNDVS" 
TYPE  STNDV(2)    ,    STNDV(l) 
TYPE   "WEEKDAY-WEEKEND  MINS" 
TYPE  RKD}(2)    ,   RMIN(l) 
TYPE   "WEEKDAY-WEEKEND  MAXS" 
TYPE  RMAX(2)    ,   RMftX(l) 
PAUSE  "PRESS  ANY  KEY  TO  CONTINUE" 
265  CONTINUE 

C***  WRITE  WEEKEND-WEEKDAY  COMPUTATIONS  TO  LINE  PRINTER 

WRITE  (12, 299) 

DO  270  I  =  1  ,  NUMSITE 

WRITE  (12,300)  SITE(I),TOWER(I),CUnNWS(I),RATEDWS(I), 
*  FRACTION(I)*100.0 

270    CONTINUE 

WRITE  (12,590)  RATED 

IF  (RCNTRL)  WRITE(12,595)  MAXRAMP 

WRITE  (12,311) 

WRITE  (12,313)  MEAN(l)  ,  MEAN(2) 

WRITE  (12,314) 

WRITE  (12,313)  STNDV(l)  ,  STNDV(2) 

WRITE  (12,315) 

WRITE  (12,313)  RHIN(l)  ,  RMIN(2) 

WRITE  (12,316) 

WRITE  (12,313)  RMAX(l)  ,  RMAX(2) 
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311 
313 
314 
315 
316 


FORMAT!"   ",///, 21X, "WEEKEND-WEEKDAY  MEANS   (ABSOLUTE) ",/) 
FORMAT!"   \14X,"  WEEKEND  — >   ",F6.2,10X,"  WEEKDAY  — >   \F6.2) 
FORMAT(///,"   ",15X, "WEEKEND-WEEKDAY  STANDARD  DEVIATIONS   (ABSOLUTE)",/) 
FORMATf///,"   ",25X,  "WEEKEND-WEEKDAY  MINIMUMS",/) 
PORMATt///,"   ",25X,  "WEEKEND-WEEKDAY  MAXIMUMS",/) 


C***  WRITE  WEEKEND-WEEKDAY  DATA  TO  DISK 

TYPE   "   " 

TYPE  "WRITING  WEEKEND-WEEKDAY  DATA  TO  DISK" 
CALL  OPENW  (1  ,  "FILENAME  FOR  MEAN  WEEK  ?  " 
CALL  OPENW  (2  ,  "FILENAME  FOR  STDV  WEEK  ?  " 
CALL  OPENW  (3  ,  "FILENAME  FOR  RMIN  WEEK  ?  " 
CALL  OFEtW   (4   ,    "FILENAME  FOR  RMAX  WEEK  ?   " 

CALL  WRITRW   (1,1,  MEAN  ,   1   ,    IERR) 

CALL  WRITRW   (2,1,  STNDV  ,   1   ,    IERR) 

CALL  WRITRW   (3,1,  RMIN  ,   1   ,    IERR) 

CALL  WRITRW   (4,1,  RMAX  ,    1    ,    IERR) 


SIZE) 
SIZE) 
SIZE) 
SIZE) 


CALL  CLOSE   (1 
CALL  CLOSE   (2 


IERR) 
IERR) 


CALL  CLOSE    (3   ,    JERR)  103 

CALL  CLOSE    (4    ,    IERR) 

450  CONTINUE 

C***  CCTPOT  PERCENT  ENERGY  LOST  IN  CCNTRCL  ALGORITHM 

IF    (.NOT.RCNTRL)    GO  TO  510 

WRITE   (12,299) 

DO  480   I  =  1   ,  NUMSITE 

WRITE   (12,300)   SITE(I),TCWER(I),COTINWS(I),RATEDWS(I), 

*  FRACTION(I)*100.0 
480            CONTINUE 

WRITE   (12,590)   RATED 

WRITE(12,595)   MAXRAMP 

WRITE (12, 500)    (WPLOSS/WPSUM)*100.0 
500  FORMAT(//,"   ",17X, "PERCENT  ENERGY  LOSS  FOR  CONTROL  =  "   ,    F5  2   "  i«l 

510  CONTINUE  '  '  ' 

C***  CCTPOT  RAMPING  DISTRIBUTION  TO  LINE  PRINTER 

WRITE  (12,299) 

DO  470  I  =  1  ,  NUMSITE 

WRITE  (12,300)  SrrE(I),TCWER(I),ajnNWS(I),RATEEWS(I), 

*  FRACTION(I)*100.0 
470    CONTINUE 

WRITE  (12,590)  RATED 

IF  (RCNTRL)  WRITE(12,595)  MAXRAMP 

WRITE  (12,351) 

DO  390  I  =  1  ,  25 

RATEL  =  FLOAT(I-1)/5.0 

RATE2  =  (FLQAT(I-l)+25.0)/5.0 

WRITE  (12,350)  RATE1  ,  FRES(I)  ,  RATE2  ,  PRES(I+25) 
390     CONTINUE 

WRITE  (12,353)  10.00  ,  PRES(51) 

WRITE  (12,299) 

DO  490  I  =  1  ,  NUMSITE 

WRITE  (12,300)  SITE(I),TCWER(I),COTI1KS(I),RATEDWS(I), 
FRACTION(I)*100.0 
490     CONTINUE 

WRITE  (12,590)  RATED 

IF  (RCNTRL)  WRITE(12,595)  MAXRAMP 

WRITE  (12,352) 

DO  395  I  -  1  ,  25 

RATE1  =  FLOAT(I-1)/5.0 
RATE2  =  (FLOAT(I-1)+25.0)/5.0 

WRITE  (12,360)  -RATE1  ,  NRES(I)  ,  -RATE2  ,  NRES(I+25) 
395     CONTINUE 

WRITE  (12,363)  -10.0  ,  NRES(Sl) 

351  FORMAT(//,"   ",20X, "POSITIVE  RAMPING  PERCENTAGE  OF  TIME",//) 
360  FORMAT("   ",4X,2(5X,  "RAMP  <  ",F6.2,"  — >  ",F6.2,"  %")) 

350  FORMATt"   ",4X,2(5X, "RAMP  >   ",F6.2,"  — >   ",F6.2,"  %"  ) 

352  FORMAT(//,"   ",20X, "NEGATIVE  RAMPING  PERCENTAGE  OF  TIME",//) 
363  FORMATt"   ",40X,"RAMP  <  ",F6.2,"  — >  ",F6.2,"  %") 

353  FORMATC  ",40X,"RAMP  >  \F6.2   "  ->  "  F6  2  "  %" 
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590  FORMAT(//,"   ",23X,  "RATED  FARM  POWER  =   ",F7.1,"  Mil") 

595  PORMAT(//,"    ",15X, "MAXIMUM  ALLOWED  POSITIVE  RAMP  =   ",F7.2,"  MW/MHJ") 

C***  OUTPUT  PERCENT  RESULTS  TO  DISK 

CALL  OPENW   (3   ,    "+RAMP  %  OUTPUT  FILENAME  >   "   ,    204   ,   SIZE) 

CALL  WRITRW    (3,1,    PRES   ,   1    ,    IERR) 
CALL  CLOSE   (3   ,    IERR) 

CALL  OPENW   (3   ,    "-RAMP  %  OUTPUT  FILENAME  >   "   ,   204   ,   SIZE) 

CALL  WRITRW    (3,1,   NRES   ,    1    ,    IERR) 
CALL  CLOSE    (3   ,    IERR) 

END 
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ABSTRACT 

The  thesis  shows  the  effects  that  wind  generated  electric 
power  has  on  conventional  electric  generation  in  Kansas  for  the 
years  1980  to  1982.  To  perform  this  analysis,  ramping  rates  in 
MW/min  of  the  utility,  wind,  and  utility  with  added  wind  genera- 
tion were  inspected  and  showed  detrimental  effects  to  utility 
operation.  However,  variations  in  cut-in  and  rated  turbine  wind 
speeds,  penetration,  wind  array  size,  and  amounts  of  upwards 
ramping  control  were  incorporated  and  showed  trends  of  decreasing 
negative  ramping  effects  due  to  injected  wind  generation.  From 
this,  it  was  found  that  a  three  site  farm  at  15%  penetration  with 
appropriate  control  parameters  behaved  with  the  same  impact  as  a 
single  site  farm  at  5%  penetration.  Also  discussed,  was  the 
correlative  relationship  between  Kansas  load  demand  and  Kansas 
wind  power  production.  For  single  site  farms,  a  direct  diurnal 
correlation  was  seen,  while  multiple  site  wind  arrays  showed 
decreasing  diurnal  correlation  for  an  increasing  number  of  com- 
bined wind   farms. 


